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CHAPTER  1 


PERFORMANCE  COMPARISON  OF  PRCPM  SCHEMES 
WITH  INTERFERING  SIGNALS 

1.1  Introduction 

Partial  Response  Continuous  Phase  Modulation  (PRCPM) 
schemes  have  recently  attracted  attention  because  of  their  improved 
performance  capabilities  over  traditionally  used  digital  modulation 
schemes.  These  capabilities  occur  from  observing  the  signal  over  more 
than  one  symbol  interval,  thereby  utilizing  the  memory  inherent  in  the 
continuous  phase  transitions. 

The  ultimate  measure  in  a  digital  communication  system  is  the 
probability  of  error.  However,  formulation  in  terms  of  the  error  proba¬ 
bility  is  feasible  only  in  simple  cases.  Often  an  intermediate  measure  is 
required  to  be  considered  and  the  error  probability  has  to  be  evaluated 
under  certain  simplifying  assumptions.  In  the  adjacent  channel 
interference  problem  the  Intermediate  measures  usually  considered  are 
the  mean-square  crosstalk  and  the  worst-case  crosstalk  [1,7].  The 
worst-case  crosstalk  and  usually  requires  extensive  simulation  to  com¬ 
pute,  while  mean-square  crosstalk  is  much  more  amenable  to 
mathematical  formulation  and  also  has  the  meaning  of  noise  variance. 

This  chapter  presents  the  signal  to  mean-square  crosstalk  as  a 
performance  measure  among  a  selected  set  of  receiver  filters  for  various 
CPM  schemes.  Also,  it  simulates  and  compares  the  error  performance 
in  the  presence  of  ACI  (with  a  different  number  of  interferers)  for  vari¬ 
ous  CPM  schemes  and  among  the  AOF  receiver  with  observation  inter¬ 
vals  of  three  and  four  bits.  In  the  next  section  we  will  introduce  the 
various  modulation  schemes  that  have  been  used  throughout  our  work 
and  emphasize  an  easy  way  to  calculate  the  power  spectra. 

1.2  Modulation  Schemes 

The  bandwidth  efficiency  of  a  modulation  scheme  is  defined  as  the 
ratio  of  the  bit  rate  transmitted  to  the  bandwidth  used.  This  is 
equivalent  to  the  number  of  transmitted  bits  per  cycle  of  channel 
bandwidth.  We  have  shown  in  Equation  1.1  and  have  repeated  here 
for  convenience 


Bandwidth  Efficiency  = 


logoM 

~BT~ 


bits/sec/Hz 


(1.1) 


that  signals  with  small  BT  products  are  the  most  bandwidth  efficient. 
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while  no  single  universal  definition  is  available  for  bandwidth. 
Amoroso  [15]  has  shown  that  power  spectral  density  is  a  key  parame¬ 
ter  in  the  definition  of  bandwidth  efficiency.  We  will  use  this  approach 
to  evaluate  numerically  the  spectrum  of  various  modulation  schemes 
as  a  measure  of  the  bandwidth  efficiency. 

The  modulation  schemes  and  the  corresponding  pulse  shapes  con¬ 
sidered  in  this  work  are  Raised  Cosine  (RC) 


g(t)  = 


The  frequency 
the  scheme  is  defined  as  LRC. 

Tamed  Frequency  Modulation  (TFM) 

6(t)  =  [6o(t-T)  +  2g„(t)  +  g„(t+T)j  (1.3) 


^(1  -  oos(|^))  ;  o  <  t  <LT 

o  :  otherwise 


(1.2) 


pulse  g(t)  is  a  raised  cosine  over  L  symbol  intervals  and 


where 


Gaussian  MSK  (GMSK) 


T 

.  .  T 

t  — 

g(t)  - 

Q 

27rBb  .  - 

Vln2  . 

-Q 

27rBb  ; 

Win  2  . 

(1.4) 


(1.5) 


where  Q(x)  is  the  Gaussian  error  function. 
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1.3  Error  Probability  Measure  in  Interfering  Environment 


The  desired  transmitted  signal  is  given  by: 


s(t,  a)  = 


2E 


—  cos(c^^t  +  0(t,  a)) 


(1.6) 


where  iJ:  =  Oq,  •  •  •  Q;n_i,  QJj,,  •  •  •  is  an  infinitely  long  sequence  of 
independent  binary  data  symbols,  taking  values  i-l.  E  is  the  symbol 
energy  and  T  is  the  symbol  time,  is  the  carrier  frequency  and  the 
information  carrying  phase  is  given  by: 


<?i(t,  a)  =  27rh  X)  q(t— iT)  (1.7) 

oo 

h  is  the  modulation  index.  The  phase  response  q(t)  is  defined  by: 

t 

q(t)  =  /  g(r)  dr 

—  OO 


where  g(r)  is  the  frequency  pulse  function. 

The  adjacent  channels  are  distributed  two  on  each  side  of  the  desired 
channel.  Due  to  the  symmetry  between  both  sides,  one  side  is  intro¬ 
duced  and  the  other  is  similar. 

The  interfering  signals  on  one  side  of  the  desired  channel  are 

^Si(t,^)  for  l<k<2.  for  l<k<2  are  the  amplitude  of  the 
k— 1 

interfering  signals  and  are  assumed  to  be  unity  among  the  k  interfer- 
ers.  The  center  frequency  of  the  first  interferer  Wj  is  Aw  rad/sec  away 
from  Wq  and  the  second  interferer  W2  =  2Wi  is  2Aw  red/sec  away  from 
Wq.  There  is  a  delay  of  r  and  2r  secs  between  the  first  and  second 
interfering  signals  and  s(t,GL)  respectively.  Also  for  1  <k<2  are  the 
random  phase  offests.  This  implies  that  the  interfering  signals  are 
independent  of  one  another  and  of  the  desired  signal.  The  interfering 
signals  can  be  represented  as 

S  si(t,4)  =  cos  fnCt-’-k)  +  ’/<t-Tk,4)  -f  f J  (1.8) 

k-l  ^  ^  k-l  ■* 


rjj  and  are  uniformly  distributed  in  the  range  [0,  T]  and  [O,  27r] 
respectively.  The  simulated  probability  of  error  for  the  whole  model  is 
given  by 


N. 


S  1-1 


iESSD  P  (^>Xj+kxii,) 

k  1  m  n 


where  k,  1,  m,  and  n  take  values  +.  1  with  equal  probability,  xjj^  for 
l<k<4  are  the  outputs  due  to  the  interfering  signals  and  are  com¬ 
puted  numerically.  W,  is  the  weight  of  the  pdf  of  the  noise  and  is 
used  to  reduce  the  computational  time  (Modified  Monte- Carlo 
method). 

The  error  performance  of  various  schemes  has  been  investigated 
with  and  without  the  effect  of  ACI  over  an  AOF  receiver  with  observa¬ 
tion  length  of  3  bits.  Figure  1.1  shows  the  theoretical  probability  of 
error  and  it  shows  that  the  MSK  scheme  performs  better  than  3RC  or 
TFM5  or  GMSK5.  The  excellent  correlation  between  the  theoretical 
results  and  the  simulated  results  in  the  absence  of  ACI  is  shown  in  the 
same  figure.  This  gives  us  confidence  about  the  accuracy  of  the  esti¬ 
mates. 

The  simulated  error  probability  with  the  effect  of  ACI  at  different 
channel  spacings  (Af)  is  illustrated  in  Figure  1.2  and  1.3  for  Af=0.8 
and  1.2  respectively.  It  can  be  seen  in  Figure  1.2  that  MSK  and 
GMSK5  with  /?T=0.2  are  severely  degraded  at  Af=0.8  till  the  point 
that  error  probability  of  10~^  was  unattainable  at  SNR=14  dB,  while 
3RC  and  TFM5  schemes  both  have  a  better  ACI  rejection  at  the  same 
channel  spacing.  At  larger  channel  spacing  (Af=1.2)  shown  in  Figure 
1.3,  MSK  has  the  best  error  performance  due  to  less  effect  of  ACI  at 
larger  channel  spacing,  3RC  and  TFM5  both  have  better  ACI  rejec¬ 
tion,  but  GMSK5  still  degrades  badly  and  error  probability  of  10"^ 
was  unattainable  at  SNR=12  dB. 

In  Figure  1.4,  we  looked  at  the  theoretical  error  t/robability  for 
CPM  schemes  at  a  larger  receiver  observation  length  (four  bits).  Both 
TFM  and  GMSK  schemes  with  (/?T=0.25)  were  truncated  over  four  bit 
periods.  It  was  found  that  for  AOF(4T),  GMSK4  outperforms  TFM4 
while  the  3RC  scheme  performs  the  best  among  all  of  them.  The  error 
probability  with  two  interferers  (k=2)  at  lower  carrier  to  interference 
ratio  (C/I=— 10)  is  shown  in  Figure  1.5  and  1.6  for  Af=0.8  and  1.2 
respectively.  One  could  note  that  the  degradation  in  SNR  increases  at 
specific  error  probability  as  the  carrier  to  interference  ratio  decreases. 

1.4  Signal  to  Mean-Square  Crosstalk  Measure 

In  this  section,  we  will  present  a  numerical  evaluation  of  the  sig¬ 
nal  to  mean-square  crosstalk  as  a  performance  measure  for  various 
modulation  schemes  and  among  different  MSK-type  receiver  filters. 
The  derivation  of  the  mean-square  crosstalk  has  been  based  on  the  fol¬ 
lowing  assumptions; 
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Figure  1.5  Error  probability  for  various  CPM  schemes  in  the  presence  of  ACI, 
(k=2),  at  lower  values  of  C/I,  Af=0.8. 
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Figure  1.6 


Error  probability  for  various  CPM  schemes  in  the  presence  of  ACI, 
(k=2),  at  lower  values  of  C/I,  Af=1.2. 


1.  The  carrier  to  interference  ratio  is  unity. 

2.  No  intersymbol  interference  other  than  the  one  caused  by  the 
modulation  scheme. 

3.  The  interfering  signal  is  stationary. 

4.  No  special  processing,  such  as  filtering  or  data  windowing,  is 
applied. 

The  mean-square  crosstalk  [8j  is  defined  as 

X  =  ?)]  (1-10) 

Where  E  is  the  expected  value.  The  signal  to  mean-square  crosstalk 
can  be  evaluated  numerically  by 


SCtR  = 


xV(d.V2)  E 


(1.11) 


where 


X  = 


^  r> 


/  cos  (B(t)  ja(t-l-r)  dt 

-kT  ''  ^ 


+ 


/  sin  lB(t)  ja(t-l-r)  dt 
-kT  ^  ' 


dr 


(1.12) 


and 

B(t)=  [Au;t-|-V<t,,^))  (1.13) 

By  using  (1.11),  we  were  able  to  compute  the  signal  to  mean-square 
crosstalk  ratio  numerically  as  a  function  of  channel  spacing.  The 
numerical  results  among  a  selected  set  of  receiver  filters  for  3RC, 
TFM(5T),  and  GMSK(5T)  schemes  [ll]  are  illustrated  in  tables  1.1 
through  1.3  respectively. 


Table  1.1.  Signal  to  mean-square  crosstalk  ratio  in  a 
CPM  system  employing  a  SRC  scheme. 


Filter 

Filter 

Channel  Spacing 

Ivind 

Length 

0.4 

0.8 

1.2 

1.6 

2.0 

MSK 

4.76 

18.10 

28.90 

35.72 

39.27 

AOF 

5T 

5.21 

21.81 

35.81 

43.10 

46.50 

MSEOF 

5T 

5.24 

22.90 

40.45 

55.45 

62.91 

AMF 

4T 

5.31 

24.09 

45.82 

66.91 

76.27 

Table  1.2.  Signal  to  mean-square  crosstalk  ratio  in  a 
CPM  system  employing  a  TFM5  scheme. 


Filter 

Filter 

Channel  Spacing 

Kind 

Length 

0.4 

0.8 

1.2 

1.6 

2.0 

AOF 

4T 

5.26 

22.37 

29.95 

32.56 

34.55 

MSEOF 

4T 

5.34 

24.63 

35.49 

47.21 

53.86 

AMF 

4T 

5.43 

25.92 

41.08 

58.37 

67.13 

Table  1.3.  Signal  to  mean-square  crosstalk  ratio  in  a 
CPM  system  employing  a  GMSK4  scheme. 


Filter 

Filter 

Channel  Spacing 

Kind 

Length 

0.4 

0.8 

1.2 

1.6 

2.0 

AOF 

5T 

5.30 

24.52 

39.02 

44.23 

47.46 

MSEOF 

5T 

5.42 

25.62 

42.08 

57.67 

65.24 

AMF 

4T 

5.76 

26.41 

47.62 

68.27 

78.13 

The  signal  to  mean-square  crosstalk  has  also  been  evaluated  by  com¬ 
puter  simulation  of  Equation  1.10.  Both  the  numerical  and  the  simu¬ 
lated  results  are  well  correlated  and  are  shown  in  Figures  1.7  to  1.9  for 
SRC,  TFM(5T),  and  GMSK(4T)  schemes  respectively. 

The  three  figures  show  the  signal  to  mean-square  crosstalk  ratio 
(SC-pR)  against  channel  spacing  among  a  selected  set  of  MSK-type 
receiver  filters  for  three  partial  response  CPM  schemes.  It  was  found 
that  AMF(4T)  and  MSEOF(5T)  receivers  have  a  substantial  improve¬ 
ment  in  SCpR  with  respect  to  AOF(4T)  and  MSK  receivers  at  all 
channel  spacings.  Figure  1.7  shows  that  the  MSK  receiver  performs 
the  worst  in  the  SCpR  measure  and  it  was  shown  in  the  previous  sec¬ 
tion  that  the  same  receiver  filter  degrades  the  most  in  the  error 
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ACI  rejection  of  MSK-type  receivers  for  a  Figure  1.8 
3RC  scheme. 


ACI  rejection  of  MSK-type  receivers  for  a 
TFM5  scheme. 
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ACI  rejection  of  MSK-type  receivers  for  a 
GMSK4  scheme. 
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probability  measure  with  interfering  signals.  While  the  AOF  receiver 
does  not  have  as  good  signal  to  crosstalk  performance  as  the  AMF 
receiver,  it  has  the  best  error  performance  among  all  other  receiver 
filters  in  the  interfering  environment.  One  last  observation  is  that 
although  the  MSEOF  receiver  is  inferior  in  error  performance  with 
interfering  signals  to  the  AOF  receiver,  it  has  a  much  better  rejection 
of  crosstalk  than  the  AOF  receiver  at  all  channel  spacings.  This  has 
been  analyzed  for  three  modulation  schemes  as  seen  in  the  last  three 
figures. 

1.5  Discussion  and  Conclusions 

In  this  chapter,  the  interferences  have  been  modeled  as  adjacent 
channel  interference  and  co-channel  interference.  The  first  model  was 
considered  with  the  effect  of  two  and  four  interferers  (k=2,4)  on  the 
probability  of  error.  It  was  found  that  the  MSK  scheme  has  optimum 
performance  in  the  absence  of  ACl  and  degrades  the  most  in  the  pres¬ 
ence  of  ACI.  It  was  also  observed  that  some  PRCPM  schemes  like 
GMSK5  degrade  as  badly  as  MSK  at  closer  channel  spacing  (Af=0.8), 
while  SRC  and  TFM5  both  have  a  very  good  rejection  of  ACI  at  all 
channel  spacings  for  AOF(3T).  One  more  important  observation;  it  is 
not  always  true  that  the  error  performance  improves  at  a  larger 
receiver  observation  length.  This  result  can  be  seen  in  Figures  1.1  and 
1.4.  So,  one  may  conclude  that  large  improvements  can  be  obtained  if 
the  detector  filter  length  is  selected  appropriately.  It  was  also  shown 
that  PRCPM  schemes  have  much  lower  sidelobes  in  relation  to  full 
response  scheme  (MSK).  These  results  have  made  PRCPM  schemes 
popular  in  bandlimited  frequency  reuse  mobile  radio  systems,  espe¬ 
cially  when  large  system  capacity  is  desired. 

It  was  found  that  the  AMF  and  MSEOF  receivers  perform  sub¬ 
stantially  better  than  the  AOF  and  MSK  receivers  for  SRC  and  per¬ 
form  better  than  AOF  for  TFM6  and  GMSK4  schemes.  This  means 
that  the  AMF  and  MSEOF  receivers  have  a  better  rejection  of  ACI  at 
all  channel  spacings  in  comparison  to  the  AOF  receiver  employing 
PRCPM  schemes.  One  concludes  that  in  order  to  choose  the  suitable 
receiver  filter  for  a  specific  application,  one  must  compromise  between 
the  error  performance  and  the  crosstalk  performance. 


I 
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CHAPTER  2 

SPACE  DIVERSITY  OF  CPM  OVER  RAYLEIGH 
FADING  CHANNELS  WITH  INTERFERING  SIGNALS 

2.1  Introduction 

In  a  mobile  radio  environment,  the  signal  transmitted  between  a 
fixed  base  station  and  a  moving  radio  unit  signal  suffers  from  fading 
and  interference.  Fading  is  caused  due  to  multipath  transmission. 
There  are  also  two  sources  of  interferences  in  this  environment.  The 
first  source  of  interference  is  the  ACI  and  the  second  source  of  interfer¬ 
ence  is  the  CCI.  These  interferences  and  fading  represent  a  major 
source  of  degradation.  Diversity  combining  is  one  of  the  techniques 
that  can  be  applied  successfully  to  combat  fading  and  reduce  interfer¬ 
ence. 

In  the  next  section,  we  will  consider  how  different  envelopes  of 
fading  signals  can  describe  the  transmitted  channel. 

2.2  CPM  over  Rayleigh  Fading  Channel 

In  this  section  an  independent  slow  Rayleigh  fading  channel  is 
assumed.  The  error  probability  for  large  signal-to-noise  ratios  is 
derived  for  CPM  over  Gaussian  and  Rayleigh  fading  channels.  Also, 
the  error  probability  is  shown  for  Rayleigh  fading  channels  with 
interfering  signals, 

2.2.1  Error  Probability  on  Rayleigh  Fading  Channels 

The  fading  signal  s(t)  is  assumed  to  follow  the  Rayleigh  density 
function  [2]  given  by 


-s»/n 


where  fl  is  the  instantaneous  signal  and  is  given  by 

fl  =  ^s^^  . 


(2.1) 


(2.2) 


The  instantaneous  SNR  can  be  written  as 


(2.3) 
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where  Nq  is  the  spectral  height  of  the  one-sided  Gaussian  noise.  By 
substituting  (2.3)  and  (2.2)  into  (2.1),  one  can  get  the  pdf  of  as 

f(i^)  =  ^e-*'/f'  ,  i/>0  (2.4) 

where  F  is  the  average  SNR  and  is  defined  by 


The  average  error  probability  of  the  modulation  scheme  with 
coherent  detection  can  be  written  as: 


m 

P(i^)  =  E  CiQ 


1-1 


(2,6) 


where  Q(  • )  is  the  error  function  associated  with  the  normal  distribu¬ 
tion  dj“  is  the  squared  Euclidean  distance  associated  with  a  signal 
corresponding  to  data  sequence  number  i  received  in  the  fixed  filter, 
and  Cj  is  the  probability  of  that  specific  signal,  m  is  the  total  number 
of  sequences.  Assuming  independent  data  symbols  with  equal  probabil¬ 
ity,  then 


'  m 


The  average  bit-error  probability  in  the  fading  case  is 


(2.7) 


P  _  /  P(j-)f(j/)di/  (2.8) 

0 

where  P(i/)  and  f(i/)  are  given  by  (2.6)  and  (2.4)  respectively.  Using 
(2.6)  and  (2.4)  in  (2.7)  one  can  obtain 


P  = 


2m 


j-i 


l+d,2r/2 


(2.9) 


From  (2.9)  one  can  evaluate  the  bit-error  probability  over  a  Rayleigh 
fading  channel  for  CPM  schemes.  Figure  2.1  shows  the  error  probabil¬ 
ity  for  SRC,  TFM4,  and  GMSK4  schemes  employing  an  AOF(4T) 
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receiver  on  Gaussian  and  Rayleigh  fading  channels. 


2.2.2  An  Easy  Way  to  Evaluate  the  Error  Performance  With 
Interfering  Signals 

In  this  section,  we  develop  an  easy  method  to  evaluate  the  error 
probability  in  the  interfering  environment.  This  method  might  not  be 
as  accurate  as  the  simulation  technique  used  previously  but  it  is  much 
faster  from  the  computational  point  of  view.  The  analysis  here  is 
based  on  stationary  white  Gaussian  noise  with  power  spectral  density 
No- 

The  algorithm  is  mainly  an  evaluation  of  the  normalized  squared 
Euclidean  distance  of  the  desired  and  interfering  signals.  For  the 
desired  signal,  the  normalized  squared  Euclidean  distance  is  calculated 
from: 


2 

T 


j  a(— t)  s(t,Q;j)  dt 
1 


/a-(t)  dt 

I 


(2.10) 


For  the  interfering  case,  the  normalized  squared  Euclidean  distance  is 
computed  numerically.  It  is  also  computed  by  computer  simulation  of 
the  filter  output  due  to  each  interfering  channel.  The  filter  output  due 
to  the  jth  interfering  channel  is 

[a(t)*cos(Wj(t-rj)  +  i/<t-rj,^)  +  fj)]  (2.11) 


The  normalized  squared  Euclidean  distance  of  the  jth  interferer  is 
defined  as  (dj^^)  and  is  given  by 


(2.12) 


The  overall  normalized  squared  Euclidean  distance  for  the  jth  sequence 
is  (dj^^^)  and  is  given  by 


1-1 


(2.13) 


J 


I 
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■where  ^  represents  the  normalized  squared  Euclidean  distance  due 
to  k  interferers.  The  error  probability  for  the  interfering  case  is 


1  m 

P(^)  =  —  E  Q 

j-i 


(2.14) 


where  Q  is  the  Gaussian  error  function  and  m  is  the  total  number  of 
sequences. 

It  was  found  that  the  resulting  error  probability  using  this  algo¬ 
rithm  is  the  same  as  the  one  obtained  by  the  simulation  technique 
used  previously  with  an  error  rate  of  less  than  5%.  This  is  seen  in  Fig¬ 
ure  2.2  and  2.3  respectively.  One  could  notice  in  both  figures  that  at 
lower  C/I,  i.e.,  when  the  desired  signal  power  is  weak  compared  to  the 
interfering  signal  power,  the  performance  degradation  is  the  highest. 
Figure  2.4  shows  the  error  performance  of  a  SRC  scheme  employing  an 
AOF(4T)  receiver  with  one  and  multiple  CCI  at  various  C/I.  It  is 
recognized  that  the  difference  in  performance  degradation  due  to  the 
interferers  is  the  most  at  low  C/I  and  it  reduces  at  higher  values  of 
C/I.  One  last  observation  is  that  the  effect  of  one  CCI  and  two  CCI  is 
almost  the  same  at  higher  values  of  C/I,  i.e.,  when  the  signal  power  is 
strong  enough  the  number  of  co-channel  interfering  signals  will  not 
have  much  effect  on  the  system  performance. 


2.2.3  Error  Performance  over  Fading  Channels  ■with  Interfer¬ 
ing  Signab 

In  this  section,  we  will  evaluate  the  error  probability  when  the 
signal  is  combined  with  fading  and  interference.  It  is  assumed  that 
fading  due  to  the  interfering  channels  is  negligible  compared  to  fading 
due  only  to  the  desired  channel.  The  average  bit  error  probability  in 
this  case  is  obtained  by  substituting  (2.4)  and  (2.12)  in  (2.8).  The 
result  is 


P  = 


2m 


m 

E 

j“i 


diov^iy2 

i-fd^^,2r/2 


(2.15) 


Figure  2.5  shows  the  error  probability  for  the  3RC  scheme  on  fading 
and  both  ACI  and  CCI  with  various  C/I.  One  may  notice,  that  the 
signal  affected  by  CCI  and  fading  has  larger  C/I  than  that  affected  by 
ACI  and  fading  even  at  closer  channel  spacings.  This  is  predictable 
because  CCI  has  more  influence  on  performance  degradation  than  ACI 
due  to  it  is  existence  within  the  same  signal  bandwidth.  Figure  2.6 
illustrates  the  same  result  for  various  CPM  schemes  under  the 


2ACI  at  Af-l  J 


Figure  2.2  Effect  of  carrier  to  interference  ratios  C/I’s 
on  error  probability  for  a  3RC  scheme 
employing  an  AOF{4T)  receiver. 


Figure  2.3  Effect  of  carrier  to  interference  ratios  C/I’s 
on  error  probability  for  a  SRC  scheme 
employing  an  AOF(4T)  receiver. 


Figure  2.4  Effect  of  single  and  multiple  co-channej 
interference  CCI,  with  various  C/I’s,  on 
error  probability  for  a  3RC  scheme  using 
an  AOF(4T). 
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Figure  2.5 


Figure  2.6 


Error  probability  P  vs.  average  SNR  F  for  a  SRC  scheme  and  vari¬ 
ous  kinds  of  interfering  signals,  employing  an  A0F(4T)  receiver. 


Error  probability  P  vs.  average  SNR  F  for  various  selected  h=l/2, 
binary  CPM  schemes  with  an  AOF(4T)  receiver. 
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influence  of  ACI  and  channel  spacing  of  0.8  the  bit  rate. 

2.3  Diversity  Approach  for  CPM  with  Interfering  Signals 

Space  diversity  is  an  effective  weapon  against  the  cochannel 
interference  encountered  in  cellular  mobile  radio  systems.  High  order 
diversity  and  strong  interference  suppression  can  be  achieved  with 
modest  hardware  complexity  by  using  diversity  combining  techniques. 

In  the  next  section,  we  will  consider  ideal  maximal  ratio  combin¬ 
ing  as  an  alternative  technique  to  combat  fading  and  reduce  interfer¬ 
ence. 

2.3.1  Error  Performance  of  CPM  for  Ideal  MRC  With 
Independent  Diversity  Branches 

For  independent  diversity  branches,  there  is  an  independent  Ray¬ 
leigh  fading  in  each  branch,  while  for  ideal  maximal  ratio  combining 
the  receiver  must  know  each  path  magnitude  and  phase  to  perform 
perfect  combining.  It  must  also  have  the  property  that  the  output 
SNR  is  the  sum  of  the  instantaneous  branch  SNRs,  i.e., 

M 

(2.16) 

k-l 


where  M  is  the  number  of  diversity  branches,  is  the  instantaneous 
SNR  at  the  kth  branch.  The  fading  signal  s(t)  at  the  kth  branch  is 
assumed  to  follow  the  Rayleigh  density  function,  i.e., 

f{s)  = s  (2.17) 

where 

^k  =  M  •  (2.18) 


The  instantaneous  SNR  at  the  kth  branch  is  given  by 


; 

■ 


A 


-f 


o 


(2.19) 


Hence,  the  p.  d.  f.  of  i/),  is 
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(2.20) 


where  Fj^  is  the  average  SNR  per  the  kth  branch  and  is  considered 
equally  among  all  branches,  i.e., 

rk  =  r.  (2.21) 


In  this  case,  the  p.  d.  f.  of  u  is  obtained  by  substituting  (2.21)  and 
(2.16)  into  5.16,  hence 


(M-1)! 


(2.22) 


for  M-branch  MRC,  assuming  independent  Rayleigh  fading  in  each 
branch.  The  average  receiver  output  SNR  is 

E{u)  =  M  r.  (2.23) 


The  bit  error  probability  for  the  MRC  diversity  case,  for  PRCPM 
with  coherent  MSK-type  reception  (3],  is  obtained  by  substituting 
(2.22)  and  (2.7)  in  (2.8)  and  is  given  by 


P  = 


2m 


E 

1-1 


1  + 


V 


"?r7r- 

H-d2r/2 


1  + 


112 


d,2r  , 


+  -^(1  +  ^)-'+ 


2!2^  '  2 
1.3.5 . (2M-3) 


. .  + 

d^r 

a. 


(2.24) 


Numerical  calculations  of  (2.24)  for  various  modulation  schemes  over 
an  AOF(4T)  receiver  are  shown  for  2  and  4  branches  in  Figure  2.7. 
One  may  notice  the  improvement  in  error  probability  versus  the  aver¬ 
age  per  branch  signal-to-noise  ratio  (F)  for  various  modulation  schemes 
by  using  more  diversity  branches. 

The  case  of  diversity  against  both  fading  and  interference  can  be 
evaluated  numerically  by  replacing  the  normalized  squared  Euclidean 


Avw«g«  SNR  In  dB  Avaitga  SNR  In  dB 

Figure  2.7  Error  probability  P  vs.  average-per-branch  Figure  2.8  Error  probability  P  vs.  average-per-branch 
SNR  r  for  various  h=l/2,  binary  CPM  SNR  P  for  various  h=l/2,  binary  CPM 

schemes  with  an  A0F(4T)  receiver  and  schemes  with  an  AOF(4T)  receiver  and 

M-branch  diversity  with  maximal-ratio  M-branch  diversity  with  maximal-ratio 

combining.  combining. 


Avwdg*  SNR  In  dB 


Figure  2.9  Error  probability  P  vs.  average-per-branch 
SNR  r  for  a  3RC  with  various  kinds  of 
interfering  signals,  M-branch  diversity 
with  maximal-ratio  combining,  employing 
an  AOF(4T)  receiver. 
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distance  (dj')  in  (2.24)  by  the  overall  normalized  squared  Euclidean 
distance  (dj^^^)  in  (2.13).  Figure  2.8  shows  the  numerical  evaluation  of 
space  diversity  with  maximal  ratio  combining  (MRC)  technique  for  the 
same  class  of  modulation  schemes  under  both  fading  and  interference 
for  2  and  4  diveislty  branches.  A  little  performance  degradation  can  be 
seen  compared  to  Figure  2.7  due  to  the  consideration  of  the  interfer¬ 
ence.  Figure  2.9  illustrates  the  performance  improvement  by  using  the 
MRC  technique  with  2  and  4  diversity  branches,  over  that  of  Figure 
2.5. 

2.3.2  Error  Performance  of  CPM  for  Ideal  SRC  with 
Independent  Diversity  Branches 

For  the  ideal  selection  ratio  combining  technique  the  combiner 
selects  the  diversity  branch  which  has  the  largest  SNR  for  bit  deci¬ 
sions.  The  same  branch  is  used  for  all  symbols  over  one  time  interval 
for  the  receiver  filter  under  the  assumption  of  slow  fading.  For  Ray¬ 
leigh  fading,  the  instantaneous  SNTi  in  the  kth  branch  has  the 
p.  d.  f.  shown  in  (2.20).  The  output  SNR  [u)  is  equal  to 

u  =  max.  (i/j,  Uo.  . . .,  ’  ^m)  {2-25) 

where  M  is  the  total  number  of  branches.  The  probability  density 
function  for  w  (2]  is 

f(i/)  =  ~  e-*'/f'  (1  (2.26) 

where  F  is  the  average  SNTl  per  branch.  The  average  receiver  output 
SNR  in  this  case  is 

=  rf;  i  (2.27) 

k-1 

which  of  course  increases  more  slowly  with  increasing  M  than  the 
corresponding  average  for  MRC  (2.23).  Sundberg  [2]  derived  the  aver¬ 
age  bit  error  probability  for  BPSK,  M  branch  diversity  with  selection 
combining.  For  CPM  the  derivation  is  similar.  The  derivation  is  con¬ 
sidered  under  the  assumption  of  coherent  detection  of  MSK-type 
receivers.  The  average  bit  error  probability  in  this  case  is  given  by 
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where  d,*,  i=l,  . .  m  is  defined  by  the  modulation  scheme  and  by  the 

receiver  filter.  Numerical  evaluation  of  (2.28)  for  the  same  set  of 
modulation  schemes  over  an  AOF(4T)  receiver  is  illustrated  for  2  and 
4  diversity  branches  in  Figure  2.10.  The  case  of  the  SRC  diversity 
technique  against  both  fading  and  interference  is  evaluated  numerically 
by  replacing  the  normalized  squared  Euclidean  distance  (dj^)  in  (2.28) 
by  (dj^^*)  in  Equation  2.13.  Figure  2.11  shows  the  error  probability  in 
the  interfering  environment  using  different  diversity  branches  with  the 
selection  ratio  combining  technique. 

In  the  preceding  section  the  analysis  has  been  based  on  the 
assumption  that  the  fading  signals  in  the  various  branches  are 
uncorrelated.  It  is  important  to  examine  the  possible  deterioration  of 
performance  of  a  diversity  system  when  the  branch  signals  are  corre¬ 
lated. 

2.3.3  Error  Performance  of  CPM  Using  MRC  with  Corre¬ 
lated  Branches  (M=2) 

In  the  derivation  of  this  part,  the  average  SNR’s  are  assumed  the 
same  in  both  branches.  Accordingly,  the  fading  envelope  will  take  the 
form  [52]; 


n(i-p)s 


|l/2 


2Vp 


^1/2 


[n(i-p) 


(2.29) 


where  p  is  the  power  correlation  coefficient  between  branches  such  that 
0<p<l.  Ii/2(-)  is  the  1/2  order  modified  Bessel  function  of  the  first 
kind. 

By  making  the  substitution  of  i/=s^/N^  in  (2.29),  one  gets 


Vp 

r{i-p)  ^ 

The  bit  error  probability  is  given  by; 


f(-)  =  ^ 


r(i-p) 


u 


11/2 


2(l-p)V^T 


‘l/2 


P  = 


1 

2m 


E 


i=l 


(1+Vp) 


(2.30) 


Figure  2.10 


Figure  2.11 
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Error  probability  P  vs.  average- per- branch  SNR  F  for  a  SRC  with 
an  AOF(4T)  receiver  and  M-branch  diversity  with  selection  com¬ 
bining. 


Error  probability  P  vs.  average-per-branch  SNR  F  for  a  SRC 
scheme  with  an  AOF(4T)  receiver  and  M-branch  diversity  with 
selection  combining. 


Figure  2.12  Effect  of  correlation  among  diversity 
branches  (M=2)  with  the  MRC  technique 
on  the  error  probability,  SRC  scheme 
employing  an  AOF(4T)  receiver. 


Figure  2.13  Effect  of  correlation  among  diversity 
branches  (M=2)  with  the  \fRC  technique 
on  the  error  probability,  SRC  scheme  with 
adjacent  channel  interference,  employing 
an  AOF(4T)  receiver. 


Figure  2.14  Effect  of  correlation  among  diversity 
branches  with  the  MRC  technique  on  the 
error  probability,  SRC  scheme  with  co¬ 
channel  interfering,  employing  an 
A0F(4T)  receiver. 
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V 


(1-V^) 


1+ 


d,2r(i-v^) 


(2.31) 


It  is  obvious  that  the  last  expression  is  valid  for  all  permissible  values 
of  p  except  at  p=Q. 

Limiting  Case  (p=0) 


By  using  the  ascending  series  of  the  modified  spherical  Bessel 
function  [14]  for  small  argument  in  the  form 

'\/ y  Vz  In+i/2(z)  -  i.3.5....(2n+l) 

■  14-  -j/2-  4- _ 

l!(2n+3)  2!(2n+3)(2n+6) 


the  bit  error  probability  can  be  developed.  One  may  notice  that  the 
argument  of  I( . )  is  zero  in  our  case,  so  the  above  expression  suits  our 
situation  perfectly.  The  average  bit  error  probability  in  this  case  is 
given  by 


P  = 


(1-P) 

2m 


E 


i-l 


2-1-3/7, 

2(1  +  1/7,)^/2 


where 


dmi-p) 


(2.33) 


(2.34) 


Figure  2.12  shows  the  effect  of  correlation,  among  diversity 
branches  with  maximal  ratio  combining  technique,  on  the  error  proba¬ 
bility  for  the  3RC  scheme.  Figures  2.13  and  2.14  illustrate  the  same 
effect  of  the  correlation  while  taking  into  consideration  the  influence  of 
the  adjacent  channel  interference  and  co-channel  interference  on  the 
error  performance  respectively. 

2.4  Conclusions 

In  this  chapter,  the  problem  of  the  interfering  channels  in  Ray¬ 
leigh  fading  environment  with  diversity  reception  for  partial  response 
continuous  phase  modulation  (PRCPM)  schemes  was  addressed. 


PRCPM  schemes  with  various  kinds  of  interfering  signals  on 
Rayleigh  fading  channels  perform  very  poorly  and  a  bit  error  probabil¬ 
ity  of  10“^  was  unattainable  at  SNR  of  20  dB.  By  introducing  the 
diversity  combining  approach  with  the  maximal  ratio  combining  tech¬ 
nique,  even  with  only  two  diversity  branches,  the  error  performance 
tremendously  improved.  By  using  diversity  with  the  selection  com¬ 
bining  technique  the  error  probability  could  be  improved  in  the  same 
manner  but  with  less  improvement  than  the  one  accomplished  by 
using  the  maximal  ratio  combining  technique.  This  means  that  the 
results  obtained  from  this  chapter  reduce  to  the  classical  results  for  the 
case  of  no  interference. 
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CHAPTER  3 


ADAPTIVE  MULTIPATH  DIVERSITY  RECEIVERS 


3.1  Introduction 

The  number  of  terminals  for  indoor  use  has  grown  considerably 
over  the  years.  An  emerging  concept  is  to  provide  wireless  communi¬ 
cation  between  these  terminals.  Indoor  wireless  communication  is 
hampered  by  a  severe  multipath  problem.  This  multipath  propagation 
manifests  as  Intersymbol  Interference  (ISI)  and  frequency  selective  fad¬ 
ing.  For  indoor  channel  applications  where  the  symbol  duration  is 
much  larger  than  the  multipath  spread,  frequency  selective  fading  is 
the  major  error  causing  mechanism.  Since  the  errors  in  symbol  detec¬ 
tion  occurs  due  to  envelope  fading  diversity  can  be  used  to  mitigate 
the  effect.  By  using  Direct  Sequence  Spread  Spectrum  (DSSS)  signal¬ 
ing  multipath  diversity  can  be  used  to  combat  fading.  The  DSSS  sig¬ 
nal  resolves  the  multipath  to  provide  uncorrelated  replicas  of  the 
transmitted  signal  at  the  receiver.  The  receiver  is  designed  to  demodu¬ 
late  each  replica  independently  to  attain  diversity  gain.  This  chapter 
is  concerned  with  the  formal  development  of  a  class  of  known  delay 
adaptive  multipath  diversity  receivers  and  the  evaluation  of  their 
p;erformance.  The  receiver  is  adaptive  in  the  sense  that  it  periodically 
measures  channel  characterizing  parameters  (path  gains  and  phase 
shifts)  and  uses  these  estimates  for  signal  detection. 

The  receiver  structure  is  derived  as  a  generalization  of  [4]  and  [5] 
for  a  frequency  selective  fading  quasi-static  channel.  This  formulation 
is  appropriate  only  for  channels  with  multipath  spread  less  than  the 
signaling  symbol  duration,  such  as  an  indoor  wireless  communication 
link.  No  restriction  is  placed  on  the  shape  of  the  multipath  delay 
power  profile  of  the  channel.  Orthogonal  and  phase  coded  signaling 
are  considered.  Finally,  an  attempt  is  made  to  obtain  closed  form 
expressions  for  bit  error  probability. 

3.2  Receiver  Structure  for  Binary  Orthogonal  Signaling 

The  received  signal,  in  the  symbol  interval  (L— l)T<t<LT, 
under  the  two  hypotheses  (Hj,  i  =  0,  1)  can  be  expressed  as  in  [13]: 

H(L)  =  H,  :  r(t)  =  ^  Ck(L)  A  p,(ti  -  kTp)  -f  z(t)  (3.1) 

k-l 


where  tj  =  t  —  (L— 1)T.  p,(t)  is  the  DSSS  code.  A  is  the  signal 
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amplitude.  z(t)  is  an  additive  zero  mean  complex  Gaussian  noise  pro¬ 
cess  with  autocorrelation  function  r)  =  S{t  —  r).  The  noise  is 
assumed  uncorrelated  with  the  message  sequence  and  the  Cj.'  s. 

(i)  Vector  Characterization  of  Received  Signal 

The  first  step  in  the  receiver  formulation  is  to  obtain  a  vector 
characterization  of  the  received  signal  over  each  signaling  interval. 
This  is  done  by  defining  the  following  set  of  orthonormal  basis  func¬ 
tions: 

k  =  1,  2 . Ne 

k=N,-t-l . 2N,  (3.2) 

k  >  2N, 


MO  = 


Ve: 


Pi(t  -  kTp) 


Po(t  -  k-N,  Tp) 

bUt) 


where  {bk(t),  k  >  2Nc}  is  an  arbitrary  orthonormal  set  of  basis  func¬ 
tions  that  are  chosen  to  make  the _ prthonormal  set 

{bij(t),  k  =  1,  2,  •  •  •  }  complete.  The  notation  k— N^,  implies  (k— 1\). 

The  definition  of  the  basis  function  is  based  on  the  ideal  assump¬ 
tion  that: 

T  fep  .^k-m)  j=j 

/  p,(t  -  kTp)  pj(t  -  mTp)  dt  =  .^.  (3.3) 

Define 


LT 

=  /  r(t)  b;(ti)  dt  (3.4) 

(L-i)T 


where  ri5(L)  is  a  projection  of  r(t)  along  b)^(t)  in  the  time  interval 
(L— 1)T  <  t  <  LT.  r]5(L)  is  a  linear  functional  of  the  received  signal 
r(t). 

Assuming  T^,  T,  ISI  can  be  neglected.  Substituting  (3.1)  and 
(3.2)  in  (3.4)  and  using  (3.3),  it  can  be  readily  shown  that  under  each 
hypothesis,  (H,,  i  =  0,1),  rij(L)  is  given  by: 
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H(L)  =  Hi  : 


rk(L)  = 


i  A\/Ej[^  Ct(L)  +  zJl 
(i— 1)  A\/E^ 

/kl 


k  =  1,  2,  ...,  Ne 
k  =  Nc  +  1,  •  ,  2Ne 

k  >  2N, 


where 


LT 

ZkL  =  /  Z(0  — ^  Pl(^l  ~ 

(L-l)T 


(3.5) 


LT 


ZkL  =  /  2(t)— ^  Po(ti  -  k-Ne  Tp)  dt 

(L-l)T  V^p 


ZkL 


LT 

J  z(t)  b'i.(t)  dt 

(L-l)T 

i 


Now,  (zkL,  k  =  1,  2,  Nj,  (zj^L  ,  k  =  Nc  +  1,  2Nc}  and 
{zjjL»  k  >  2N<,}  are  noise  terms,  which  are  uncorrelated  for  all  k,  L  and 
are  identically  distributed  zero  mean  complex  Gaussian  random  vari¬ 
ables  with  variance  N,,.  | 

Due  to  the  independence  of  noise  terms,  z^l^  ZkL  Zj^l,  and 
since  the  noise  terms  are  uncorrelated  with  the  Cj^,  it  is  seen  from  (3.5) 
that  {rk(L),  k  =  1,  2,  ...,  2Nc}  form  a  set  of  sufficient  statistics.  This 
set  of  sufficient  statistics  is  defined  vectorially  as 
R(L)  =  lr,(L)  r2(L)  •  ■  •  r,N,(L)lf. 

(ii)  Definition  of  Memory  Information 

It  is  proposed  to  use  a  memory  to  store  the  received  signal  over 
Mj,  past  symbol  durations.  This  means  that  at  time  t  =  LT  all  the 
information  available  to  the  receiver,  on  the  complex  path  gains  C)^,  is 
contained  in  r(t)  received  over  (L  —  Mb)T  <  t  <  (L— 1)T.  As  an  alter¬ 
native,  it  is  possible  to  store  a  linear  functional  of  r(t)  such  as  ri^(L). 
rij(L)  preserves  the  Gaussian  nature  of  the  memory  information  and  all 
the  information  on  C^.  From  (3.5)  it  is  clear  that  one  needs  to  store 
only  the  set  of  sufficient  statistics  {ri^(j),  k  =  1,  2,  ...,  2Nj,}  over 


j 
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(L  —  Mb)  <  j  <  (L~l)-  It  is  further  observed  that  each  rj^Q)  contains 
an  unknown  modulation  component,  viz.,  which  hypothesis  Hj  or  is 
true.  This  unknown  component  is  removed  from  the  memory  informa¬ 
tion,  by  redefining  the  basic  memory  variable  as: 


(i) 


VkO)  ^ 


rkO)  + 

2 


(3.6) 


(ii)  Vk(j)  4  dj  rk(j)  +  (1  -  dj)  rj^^  (j) 


(3.7) 


where  dj  =  i  when  hypothesis  H,  i  =  0,1,  is  true  in  the  signaling 
interval.  This  is  a  decision  feedback  scheme,  where  dj  corresponds  to 
past  bit  decisions.  In  (4.5)  the  division  by  2  is  done  so  that  the  Vi.(j) 
under  both  schemes  has  the  same  noise  variance  of  N^. 

Now,  the  set  of  variables  {Vjj(j),  k  =  1,  2,  ...,  N^,}  collected  over 
the  interval  (L  —  Mb)  <  j  <  (L“l)  constitutes  the  memory  information 
that  will  be  stored.  Vij(j)  are  complex  Gaussian  random  variables. 

Since  the  r(.(j)  are  uncorrelated,  the  corresponding  V)j(j)  are  also 
uncorrelated.  Hence,  the  memory  information  can  be  partitioned  into 
Nj;  independent  blocks,  each  of  length  Mb-  The  k'’^  block  of  memory 
information  is  designated  vectorially  as: 


■"kL  -  |Vk(L  -  Mb  +  1)  •  •  ■  Vj(L-1)|T  (3.8) 

The  total  memory  information  is  designated  as  where 

II"  =  [miL  m.^L  •  •  • 

Now  all  the  information  on  the  complex  channel  coefficients, 
represented  vectorially  as  C{L)  =  [C'i(L)  C^iL)  ■  ■  ■  Cj^^J^L)]  is  con¬ 
tained  in  I^.  Thus,  the  information  C(L)  given  1^*  can  be  expressed 
by  the  conditional  density: 


P{C(i)  IlD  =  n  nCM  |m,i)  (3.9) 

1 


This  expression  follow's  from  the  uncorrelated  scattering  assumption 
for  the  Ck  and  the  independence  of  memory  blocks.  Since,  is  a 
complex  Gaussian  random  vector  and  it  can  be  shown,  [l],  [24],  that: 
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P(Ck{L)|m,,) 


^  \C,{L)  -  CiL)\\ 

exp( - - ) 


(3.10) 


where, 


Ck(L)  =  E[Ct(L)  |mkJ 


(3.11) 


C)jl  is  the  MMSE  estimate  of  C)j(L)  given  the  information  vector 
and 

2<Tl\  =  E|  |C,(L)  -  Ck(L)  P  I  (3.12) 


is  the  mean  squared  error  in  estimation,  conditioned  upon  Ci^(L) 

is  a  zero  mean  complex  Gaussian  random  variable  with  variance 
2q^1^  =  Tp  <i>c(kTp:LT);  where  is  the  delay  power  profile  of  the 

channel.  The  estimate  Cj;(L)  is  also  a  zero  mean  complex  Gaussian 
random  variable  with  variance  denoted  by  2X;^^..  The  error  in  estima¬ 
tion  is  given  by: 

ek(L)  =  Ck(L)  -  Ck(L)  (3.13) 


eij(L)  can  be  easily  shown  to  be  a  zero  mean  complex  Gaussian  random 
variable  with  variance  given  by: 

=  2aL-k  -  2Xj^k  (3.14) 


(iii)  Likelihood  Ratio  Test  (LRT) 

In  formulating  the  LRT  a  Bayes  strategy  is  considered  and  the 
performance  criterion  used  is  minimum  probability  of  error.  That  is, 
if  Cl  denotes  the  event  of  a  decision  error  in  the  L*'^  symbol  interval,  it 
is  required  to  minimize: 


P(eL)  =  IlDl 


(3.15) 


Minimizing  (3.14)  is  equivalent  to  minimizing  P(eL  |Il").  The 


corresponding  LRT  is  obtained  as  a  simple  extension  from  [23].  [4]  and 
can  be  written  as: 


A(L  |lf )  - 


Ai(L  Ilf) 
A<,(L  Ilf) 


H(L)  -  H, 
HlLf-  Ho 


(3.16) 


w’here 

A,(L  |If )  -  /  P(R(L)  |H(L)  =  H„  If,  C(L)) 


•  P(C(L)  lH(L)  =  H„  If)  dC(L),  i  =  0,1  (3.17) 


Oli^ervntions; 

1.  C(L)  is  independent  of  the  hypothesis.  Further,  from  (3.6)-(3.9). 
it  is  clear  that,  P(C(L)  |H(L)  =  H,  If)  =  P(C(L)  |lf ). 

2.  The  hypotheses  are  independent  of  the  memory  information. 
This  follows  from  (3.6)  and  (3.9). 

3.  R(L)  is  independent  of  If.  This  is  easily  seen  from  (3.5).  For 
R(L)  given  C(L)  and  the  hypothesis  H(L),  the  only  uncertainty 
in  R(L)  corresponds  to  the  noise  terms  of  each  rjj(L).  These  noise 
terms  are  and  z[?j,  which  are  uncorrelated  for  all  k  and  j. 
Hence,  R(L)  is  independent  of  If. 

Using  these  observations  (3.16)  can  be  expressed  as: 


A,(L  |lf)  =  /  P(R(i)  •  C(L)) 

•  P{C(L)  I  If)  dC{L),  i  =  0,  1 


Now  from  (3.4)  and  the  definition  of  R(L)  it  is  easily  seen  that: 


P(R(L)  1  H(L)  =  H„  C(L))  = 


( — - — )'^'  n  exp[ 
^  ''  k-i  “ 


|r,(L)  -  i  aV^  Ck(L)  P 


(3.18) 
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2Ne 

n 

k-Nc+l 


exp[- 


N„ 


(3.19) 


where  i  =  0,  1.  By  using  a  binary  symmetric  hypothesis  formulation 
for  the  problem  it  was  not  necessary  to  divide  numerator  and  denomi¬ 
nator  of  the  LRT  by  the  hypothesis  under  which  only  noise  is  received, 
as  in  [4]. 

Substituting  (3.18)  and  (3.19)  in  (3.16)  the  LRT  can  be  simplified 
to  yield  the  follov/ing  receiver  structure: 


A-Ep 


No  +  A- 


Ep  2ack 


‘k^L)  |-  -f 


No 


No 


+  A-’Ep  2ct4 


2Re{v^{L)  Ck*(L)} 


H{L)  -  H, 
H(L)  -  Ho 


A% 


N. 


+  A^Ep  2ctI^ 


l>f(L)  1-  + 


No 


Np  -f  A^Ep  2cr4 


2ReK(L)  Ck*(L)} 


(3.20) 


where  r^^L)  =  rk(L)  and  r,f(L)  =  r^(L) 

Each  side  of  the  LRT,  described  by  (3.20),  is  illustrated  in  Figure 
3.1.  On  each  diversity  branch  there  is  a  dual  combination  of  coherent 
and  non-coherent  processing,  the  choice  of  one  or  the  other  is  deter¬ 
mined  by  the  factor  A^EpN^f^  2(7^^,  which  is  determined  by  the  aver¬ 
age  SNR  on  each  branch.  This  receiver  is  considered  to  be  partially 
coherent. 

Case  2  Coherent  Multipath  Diversity  Receiver 

If  the  estimation  process  yields  perfect  estimates  then  2crj^i^  =  0  and 
Cij(L)  =  Ci.(L).  The  corresponding  receiver  structure  can  be  deduced 
from  (3.19)  to  be: 
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Nc  H(L)-H,  Nc 

s  ReW(L)  C,*(L)}  E  ReK(L)  C^*{L)}  (3.21) 


This  is  similar  to  the  ideal  coherent  RAKE  correlator  [16],  [18]. 

Ca.9p  3  Non-Coherent  Multipath  Diversity  Receiver 

Consider  the  case  when  the  estimation  process  yields  very  poor  esti¬ 
mates.  Since  2(7^]^  =  2aLk  —  the  maximum  value  of  2aLi^  occurs 

when  2X]^jj  =  0.  This  implies  no  estimation  or  Cjj(L)  =  0.  Substituting 
Ck(L)  =  0  in  (3.20)  yields; 


No  H(L)-H,  No 

S  4  1%,  ?  „  E  4  |r«L)P  (3.22) 

k-l  »(L)-Hok=l 


where 


A^Ep  Tp  ■  'i».(kTp;LT)  _  X|,(LT) 

"  A^Ep  Tp  •  i’pCkTpiLT)  +  Np  Xk(LT)  +  Np 
and  from  (3.3), 

T 

riKL)  =  /  r(t)  p,(t  -  kTp)  dT  i  =  0,  1 
0 


Xjj(LT)  is  simply  the  eigenvalue  of  the  k''^  path  in  the  signaling 
interval. 

Since  DSSS  codes  are  used,  the  receiver  structure  of  (3.20),  (3.21) 
and  (3.22)  can  be  modified  for  easier  implementation.  That  is,  in  Fig¬ 
ure  3.1  the  branch  multipliers  can  be  replaced  by  switches  and  the 
integrators  replaced  by  summers.  The  summer-square  combination 
may  also  be  realized  using  an  envelope  detector,  [12]. 


3.3  Elstimation  Schemes 

From  the  receiver  structure  of  (3.20)  it  is  clear  that  for  coherent 
processing  on  each  branch  it  is  necessary  to  estimate  Cj^.  In  this  sec¬ 
tion  decision  directed  estimation  (DDE)  and  non-decision  directed  esti¬ 
mation  (NDDE)  are  considered.  The  estimator  is  formulated  based  on 
the  discrete  \Viener-Hopf  filter.  Since  the  memory  information,  used 


to  estimate  the  channel  coefficients  C^,  can  be  partitioned  into 
independent  blocks,  one  for  each  channel  coefficient,  the  estimator  need 
be  formulated  only  for  a  single  path.  All  other  branches  of  the 
receiver  will  employ  identical  estimators. 

From  (3.5)-(3.7)  it  is  clear  that  the  estimation  problem  comprises 
estimating  a  Gaussian  parameter  in  AWGN.  The  memory  information 
from  which  the  estimate  is  derived  is  also  Gaussian.  For  this  Gaussian 
estimation  problem  the  optimal  MMSE  estimator  is  linear  and  of  the 
form,  [l]: 


Ck(L)  = 


L-l 


E  w,jV,(j) 

i-L-Mb 


(3.25) 


where  Vj.(j)  is  determined  from  (3.G)  and  (3.7). 

The  weights  {W|.j  ,  j  =  (L  —  M^),  •  •  *  ,  (L— 1)}  are  chosen  to  make 
Cjj(L)  a  MMSE  estimate.  The  optimal  weights  can  be  shown,  from  the 
orthogonality  principle,  to  be  a  solution  of  M^,  discrete  Wiener-Hopf 
equations,  [1],  given  by: 


^c.(L-i)-  s'  WkilAVi;  +  hN„  i(j-i)|  , 

j-L-Mb 


i  =  (L-Mb),  •  •  •  ,  (L-2),  (L-l)  (3.24a) 


where. 


«c.(j-l)  ^  E!Ck(j)  Ck*(i)|  =  Tp  ■f>p(kTp;j-iT)  (3.24b) 


and  h=l  for  DDE  and  h=2  for  NDDE. 

The  mean  square  error  in  estimation  can  be  determined,  as  in  [14],  to 
be: 


2<^Lk=«‘c,(0)-AVE;  E  Wkj«Sc,(L-j)  (3.25) 

J-L-Mb 


For  a  quasi-static  channel  the  solution  to  (3.24a)  can  be  obtained,  for 
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any  M^,,  using  the  Sherman  Morrison  Inverse  matrix  identity,  [11]. 
The  solution  for  the  optimal  weights  can  be  shown  to  be: 


Wki  =  W'ko  = 


=  =  Wk  = 


aV^  <^Ck(0) 


A2Ep.^c,(0)  Mb  -b  hN, 


(3.26) 


From  (3.25)  and  (3.26)  the  mean  square  error  is: 


<^c.(0)  hNo 

A^Ep  <f>cj,0)  Mb  +  hN, 


(3.27) 


When  A'Ep  <ji>cj.(0)  Mb  ^  hN^,  which  occurs  when  Mb  is  large  and  the 
average  SNR  per  path  is  high,  then  (3.26)  and  (3.27)  reduce  to: 


Wk 


1 

AV^Mb 


and 


1 

A^Ep 


hN, 

Mb 


(3.28) 


(3.29) 


Hence,  when  the  average  SNTl  per  path  is  high  (weak  fade  or  large  Mb) 
(3.28)  can  be  used  in  (3.23)  to  give; 


Ck(L)  = 


1 

AVl^Mb 


E  Vk(j) 


j-L-Mb 


(3.30) 


Thus,  from  (3.30)  it  is  seen  that  the  estimate  is  derived  by  averaging 
the  memory  information  over  Mb  symbol  intervals.  This  process  is 
equivalent  to  filtering  the  memory  information  using  a  band-pass  filter 
with  bandwidth  greater  than  the  fade  rate. 

With  DDE  if  a  bit  decision  is  incorrect  it  will  propagate 
throughout  the  memory  length  causing  a  string  of  estimation  errors 
and  hence  further  bit  decision  errors.  The  effect  may  lead  to  a  rmi- 
away  condition,  [17],  and  can  be  ameliorated  by  using  a  short  memory 
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or  exponential  weights. 

The  function  of  exponential  weights,  in  the  sample  mean  estima¬ 
tor,  is  to  emphasize  the  recent  bit  decisions  with  respect  to  the  past  bit 
decisions.  One  form  of  these  weights  is: 


W'k,  =  2 


Mb-1 


Mb+1 


W, 


kj 


(3.31) 


The  use  of  exponential  weights  given  by  (3.31)  is  suggested  for  a  prac¬ 
tical  sample  mean  estimator  using  DDE.  It  can  be  shown  that  by 
using  these  exponential  weights,  there  will  be  a  small  reduction  in  the 
amplitude  of  the  estimated  path  coefficient  and  hence  a  small 
reduction  in  the  estimated  SNR  per  path. 

3.4  Performance  of  Partially  Coherent  Receiver 

3.4.1  Uniform  Delay  Power  Profile 

The  receiver  structure  of  (3.20^  can  be  modified  by  adding 
N2(a2Ep  (2o-^k)2  -1-  N^2a4)-*  |Ck(L)  P,  k  =  1,  2,  ...,  N,  to  both  sides 
of  the  LRT.  The  resulting  structure  can  be  rewritten  as: 


S  4  I  IXk(L)  P  -  |Y,(L)  r-]  0 


where 


4  = 


No 


A^Ep  2ct4  -H  N, 


No 


[A-Ep  2C7C,  J 


(3.32a) 


(3.32b) 


Xk(L)  = 


A2E„  2a?, 


Lk 


No 


t^{L)  +  Ck(L) 


(3.32c) 


Yk(L)  = 


A^Ep  2<7l\ 
No 


rk"(L)  +  Ck(L) 


(3.32d) 
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For  a  uniform  delay  power  profile  all  channels  are  identical  and 
2Q]Lk  =  k=l,  2,  Nj.  Further,  from  (3.27)  or  (3.29), 

2olj.  =  20’^.  Recalling  that  the  order  of  diversity  N(j  is  usually  less 
than  equal  to  N^.,  (3.32a)  can  be  rewritten  as: 

,  Nd  H{L)-Hi 

D(L)  &  E  IXk(L)  -  |Yk{L)  P  0  (3.33) 


The  probability  of  error  is  given  by: 

P(e)  =  P(D(L)  <  0  I  H(L)  =  H,)  P(H(L)  =  H,) 

+  P(D(L)  >  0  I  H(L)  =  Ho)  P(H(L)  =  H„)  (3.34) 


For  a  binary  symmetric  hypothesis  problem  with  equally  likely 
hypothesis,  (3.34)  reduces  to: 

P(e)  =  P(D(L)  <  0  I  H(L)  =  Hj)  (3.35) 


To  determine  P(e),  first  the  conditional  error  probability  P(e  [II")  is 
determined.  This  is  of  the  form: 

P(e  Ilf)  =  P(D(L)  <  0  I  H(L)  -  H„  If)  (3.36) 


To  evaluate  (3.37)  one  needs  to  know  the  statistics  of 
{Xi,(L)  I  H(L)  =  Hi,  If’}  and  {Yk(L)  |  H(L)  =  U,,  I£"}.  From  (3,5), 
(3.32c)  and  (3.32d)  it  can  be  shown  that  {X};(L)  |  H(L)  =  Hi,  1“}, 
k  —  1,  2,  ...,  Ng  are  i.i.d.  complex  Gaussian  random  variables  with 
mean  and  variance  (t^,  where: 


=  (1  +  aVE^Tl)  Cic(L) 

(3.37a) 

+  No) 

(3.37b) 

where 
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iL  = 


A%  2al 

N. 


Similarly,  {Yk(L)  ]  H(L)  =  Hj,  I™},  k  =  1,  2,  Nj,  are  i.i.d.  complex 
Gaussian  random  variables  with  mean  variance  Oy,  given  by: 


A'Vk  =  Ck(L) 


(3.37c) 


(3.37d) 


Also,  the  covariance,  between  the  two  conditional  variables,  can  be 
shown  to  be: 


^XkYk  ~  k  —  1,  2,  ....  Njj 


(3.37e) 


The  general  quadratic  form  of  a  decision  variable  is  given  by: 


D  =  E  lA  |Xt  P  +  B  |Yk  p  +  C  XtYk*  +  C*Xk*Yk|  (3.38) 

k-1 


where  A,  B,  C  are  constants,  {Xjj,  Yj^}  are  a  pair  of  correlated  complex 
Gaussian  random  variables  and  the  N  pairs  {Xj^,  Yj^}  are  mutually  sta¬ 
tistically  independent  and  identically  distributed. 

For  the  general  quadratic  form  of  (3.38)  a  probability  of  error 
expression  has  been  derived  in  Appendix  4B  of  [18].  From  (3.33)  and 
(3.38)  it  is  clear  that  {D(L)  |H(L)  =  Hj,  is  a  quadratic  form  in 
complex  Gaussian  random  variables  with  A  =  1,  B  =  — 1  and  C  =  0. 
Applying  the  results  of  Appendix  4B  of  [17],  the  conditional  probabil¬ 
ity  of  error  for  (3.36)  can  be  expressed  as: 

P(%  Hf )  =  QlVk?7,  Vi^l 


+  [A(Nd)  -  1]  Io(^  Xl)  exp  j^- 


ka  +  k. 
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(3.39) 


Nd 


+  S  B(n,  Nd)  I„(^  xl)  exp 


n»l 


kg  +  k4 


where 


k.Q  = 


No(2  +  7l)'7l 


,  k4  =  (A^Epn'L  +  1)^  kg 


0  =  Vi^  ,  X  =  ^^4  =  S  |Ck(L)  P 

^  2Xl  k-1 


Nd-1  2Nd-l  (1+7,)'' 


B(n,Nd) 


Nd-n-l  2Nd— 1 

S  (  k  ) 

k-0 


(1  -h  (l+Ayi^-yL)" 

(2  + 


(1  +  7l)'^^  '  (1+aVe;7l)~" 

(2  + 


Further,  !„(•)  is  the  order  modified  Bessel  function  of  the  first  kind 
and  Q(a,b)  is  the  Q-function  as  defined  in  [22]. 

The  probability  of  error  can  be  determined  by  averaging  (3.39) 
with  respect  to  the  p.d.f  of  Xl-  Given  the  memory  information  xl 
is  known.  Hence,  the  probability  of  error  is  given  by: 

OO 

P(e)  =  /  P(e  Ilf)  P(xi,)  dxi,  (3.40) 

0 


For  a  uniform  delay  power  profile  Cjj(L)  are  a  set  of  i.i.d.  zero  mean 
complex  Gaussian  random  variables  with  variance 
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2XL“k  =  2X2,  k  =  1,  2,  N,. 

Therefore,  Xl  has  a  gamma  p.d.f.  given  by: 


P(xl)  = 


N' -1 

-V'  u 


(Nd-1)!  (2\£f‘‘ 


Xl 

2^l  J 


Xl  >0 


(3.41) 


Using  (3.39)  and  (3.41)  in  (3.40)  and  performing  the  integration  yields: 


1 


X  •  2X^2 


/  \ 

-  Nh  { 

1-^ 

2F1 

x^ . 

1-N. 


.  1 4- •  1- ^ 


X- 


N<,-1 

+  E 

n— J 


Nd+n— 1 

B(n,Nd)  (  n  ) 


A 

2X 


2F1 


n+l-Nd  n-Nd  /?2 

- 0 - ’  1  +  — ^  ;  n+l; 

2  2  ^  j 


(3.42) 


where  2F i(a;  b;  c;  z)  is  the  Gauss  Hypergeometric  function  and 

00 

I(Nd)  =  /  Q(Vk^.  Vk^)  Xl exp 
0 


It  is  difficult  to  obtain  a  closed  form  expression  for  (3.43)  except  in  the 
form  of  an  infinite  series.  The  infinite  series  form  is  given  by 


XL 

2Xt2 


dxi 


(3.43) 


I(Nd)  =  E 

j-0 


/  \ 
ks 

1/2 

1 

+ 

f  \ 

1-^ 

k4 

V  ^ 

2X 

V  / 

j! 

x^J 

K 
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oFi 


j  +  1  -  N 


/?- 

;  j  +  1;  ^ 
X-  J 


(3.44) 


Using  Pad’s  algonthm,  [15],  to  compute  Q(a,b),  (3.54)  can  be  con¬ 
veniently  evaluated  by  numerical  integration  using  the  NAG  (Numeri¬ 
cal  Algorithm  Group)  routine  DOlAMF. 

3.4.2  Non-Uniform  Delay  Power  Profile 


For  a  non-uniform  delay  power  profile  the  average  SNTl  on  each 
path  is  different.  In  this  section,  a  probability  of  error  expression  is 
developed  for  a  quasi-static  channel.  For  a  quasi-static  channel  and 
the  estimator  of  (3.30)  the  estimation  error  variance  is  given  by,  (3.20): 


oack  =  2ai:  = 


hN. 


A-E_Mb 


(3.45) 


Following  the  procedure  outlined  for  a  uniform  delay  power  profile, 
the  conditional  error  probability  P(e  j  I^)  is  once  again  given  by 
(3.39).  However,  for  a  non-uniform  delay  power  profile  the  p.d.f  of  Xl 
given  by: 


Nd 


TT, 


Lk 


P(xl)  =  E 

k-I  --^Lk 


2XA 


exp 


2X4  J 


Xl>0 


(3.46) 


where 


Nd 

^Lk  =  n 
1-1 
i»‘k 


2Xr2 


Lk 


2X4  -  2X/ 


Substituting  (3.39)  and  (3.46)  in  (3.40)  and  solving  the  integral,  yields: 


P(e)  = 


Nd 

E 


k-l 


^Lk 

2X4 


A(Nd)  -  1 
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Nd-l 

+  E 


B(n,N,)  3- 

Vx=  -  0-  \\  +  Vx=  -  /?t 


(3.-1C) 


where,  x  =  (ka  +  ^4)/-  +  (^Lk)  *•  kg,  k4,  A(Nd),  B(n,Nd),  /?  are  as 
defined  in  (3.39)  and 


Ii(Nd)  =  /  Q(Vk3XL’  Vi^)  exp 


dxL 


j 

(3.-47) 


A  useful  computational  form  of  (3.-47)  can  be  obtained  by  expressing  it 
in  an  infinite  series  form  given  by: 


Ii(Nd)  -  S 

]-0 


/o 


0> 


-  0-\x  +  Vx~  - 


(3.48) 


This  is  a  convergent  series  which  can  be  evaluated  accurately  and 
rapidly,  using  the  NAG  librarj'  routine  C06B.\F,  using  a  50  term 
approximation. 

All  the  probability  of  error  expressions  derived  in  this  section  hold 
only  for  >  2Xlj..  This  is  reasonable  since  if  the  reverse  were  true 
it  would  imply  that  the  error  variance  were  negative. 

3.5  Approximations  and  Bounds 

This  section  discusses  upper  and  lower  bounds  on  the  probability 
of  error,  as  a  function  of  the  knowledge  of  channel  states.  This 
knowledge  of  channel  states  is  expressed  in  the  definition: 

E1C|,(L)  Ci*(L)) 

Ac  =  - : - (S.cig) 

lVar{0^{L)}Var{Ck(L)}|''  = 


where  k  =  1,  2,  ...,  N^j.  expresses  the  degree  of  correlation  between 
Cj.(L)  and  C|^(L).  Further,  (3.49)  can  be  simplified  and  expressed  as: 
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V 


(3.50) 


3.5.1  Lower  Bound 

When  Pk  ~  L  '*•'  follows  from  (3.14)  and  (3.50)  that  2cr£^^  =  0. 
This  implies  that  Ci^(L)  =  Cij(L),  indicating  perfect  estimation.  When 
this  condition  is  satisfied  for  all  k  branches,  then  as  in  (3.21)  perfectly 
coherent  detection  is  possible  and  the  corresponding  error  probability 
will  be  a  lower  bound. 

From  the  estimator  structure  of  (4.40),  2(t^^ — ►O  for  all  k 
branches  when  the  average  SNR  per  paih  is  high  (M^,  fixed,  No<^l). 
Under  this  high  SNR  condition  the  following  approximations,  with 
regard  to  (3.42)  hold; 

High  SNR  ApproviTri?if ionsr 

As  2(tI^  —►0,  7l  — ^  0.  B(n.  Nj)  — ►  0,  x  »  L  {—f"  -*  0,  and 
2Xl  — *•  2q'£. 

Using  these  approximations,  (3.42)  can  be  reduced  to; 

P(e)  ~ 


OO 

o  I  V^)  XL exp 

(Nj-l)!  2ql  0 


(3.51) 


J 


Now  as  “*■  0,  k3  1,  k^  1  and  k^  — *■  k3.  Also,  k^  »  k^— k3. 
Under  these  conditions  the  following  approximation  is  valid,  [22]; 


QlVkaXL-  VMl)  -  Y  erfc 


( Vi^  -  Vkl) 
\% 


ji 


j 


—  erfc 

9 


A‘EpXL 


«/2 


2  N„ 


(3.52) 


1 
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Substituting  (3.52)  in  (3.51)  and  integrating  by  parts,  yields: 


P(e)  ^  ^ 


Nd-i  2k 

1-/^1  E  (k) 

k*«0 


1  -  /V  f 


(3.53) 


where 


= 


% 


%  +  2 


l‘/2 


Ic  = 


A‘Ep  2al 

N. 


Similarly,  for  a  non-uniforni  delay  power  profile  the  lower  bound  on 
the  probability  of  error  can  be  obtained  by  making  the  same  high  SNR 
approximations.  These  approximations  allow  (3.46)  to  be  expressed  as: 


Nd 

P(e)  ~  x; 
k-1 


^Lk 

2Xi?k 


00 

/  —  erfc 
0 


A^EpXL 


2N„ 


lV2 


exp 


2X4  ) 


dxT 


(3.54) 


Integrating  (3.54)  by  parts  and  simplifying  yields: 


1 

P(e)  ^  ^  E 

k=i 


A^Ep  20^ 

1/2' 

A-Ep  2ql\  -h  2N, 

(3.55) 


Equations  (3.52)  and  (3.55)  only  hold  for  high  SNR  and  are  similar  to 
the  probability  of  error  expressions  obtained  for  a  perfectly  coherent 
RAKE  receiver,  [18]. 

3.5.2  Upper  Bound 

When  Pij  — ►  0,  the  estimate  Cj^(L)  is  totally  decorrelated  from 
Ci5^(L)  and  no  useful  estimate  can  be  obtained.  This  condition  typi¬ 
cally  occurs  when  the  channel  fades  rapidly,  as  in  a  mobile  radio  chan¬ 
nel.  This  implies  that  the  estimator  cannot  track  the  channel 
coefficients.  With  regard  to  the  estimator  of  (3.40),  this  means  that 
the  fade  rate  is  much  higher  than  the  estimator  filter  bandwidth  and 
hence  no  sample  mean  estimate  can  be  obtained.  Under  this  condition, 
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non-coherent  detection  is  done.  Although  — ►  0  can  occur  for  low 
and  high  SNR  our  interest  is  in  the  high  SNR  case  as  this  represents 
the  worst  case  and  hence  the  corresponding  error  probability  will  be 
upper  bound  on  performance. 

Equation  (3.42)  can  be  rewritten,  using  (3.44),  as: 


P(e)  = 


lA(Nd)  -  1|  ,F, 


1-Nd  ,  Nd., 

2  ’  '  2  ’ 


Nd-1  N(j  +  n— 1 


n— 1 


B(n,  Nd)  + 


2^1 


n+1  -  N’d  n  -  N’d 

- ^ .  1  +  — n  +  1;  ^ 

-  -  X-  J 


Nd  +  n 


+  E  ( 


n 


1 

) 


2F1 


n+1  -  Nj  n  -  Nj  0'. 

— 2 — ’  *  ~ — ’  ’  ly 


(3.56) 


For  high  SNTi,  corresponding  to  P(e)  >  10~^,  if  0  for 

k=l,  2 . N'd,  then  2X^1^  — ♦  0  and  =  2(tI^  for  k=l,2,...Nd. 

Further,  X  »  1,  X  ^  0  and  k3  ^  k^. 

Applying  these  high  SNR  approximations  to  (3.56)  yields: 


Nd-i  Nd  +  n  - 


n— 1 


1 

) 


B(n,Nd)  + 


(3.57) 


Equation  (3.57)  is  found  to  be  an  excellent  approximation  to  the  error 
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probability  of  a  Square-Law  receiver.  This  is  an  upper  bound  which 
also  holds  for  a  non-uniform  profile  due  to  the  following  reason.  As 
Pk  0,  2XLk  —*■  0  and  hence  2QLk  =  2cr|Lk  hNQ(A"EpMb)“\  for  all  k 
branches.  Hence,  this  case  reduces  to  that  of  a  uniform  delay  power 
profile,  since  all  paths  now  exhibit  equal  average  SNR. 

This  result  provides  an  interesting  interpretation.  When  the  esti¬ 
mate  is  completely  decorrelated  from  the  actual  channel  coefficient 
there  is  no  knowledge  of  the  path  gain.  Therefore,  the  receiver 
assumes  identical  paths,  even  though  the  profile  may  be  non-uniform, 
and  performs  Square-Law  detection.  Since  the  Square-Law  receiver 
provides  an  upper  bound  on  performance,  this  action  of  the  receiver 
may  be  viewed  as  a  minimax  operation. 

3.6  Performance  Evaluation  for  Orthogonal  Signaling 

The  performance  of  the  partially  coherent  receiver  is  evaluated  for 
an  illustrative  indoor  channel  exhibiting  an  RMS  multipath  spread  of 
250  ns.  A  data  rate  of  32  kbps  and  a  code  length  of  255  chips  were 
assumed.  Figures  3.2  and  3.3  illustrate  the  improvement  in  perfor¬ 
mance  with  higher  orders  of  diversity.  Figures  3.2  and  3.3  also  indi¬ 
cate  the  superiority  of  DDE  over  NDDE.  The  improvement  for  DDE 
is  due  to  3  dB  lower  noise  power  in  its  phase  estimation  process.  It 
has  been  assumed  that  the  DDE  scheme  is  ideal.  However,  in  practice, 
DDE  cannot  be  used  at  low  SNR  due  to  decision  errors.  NDDE  can  be 
used  at  low  SNR  and  the  receiver  can  switch  to  DDE  at  a  higher  SNTt 
or  when  there  is  a  low  error  rate.  The  switch  can  be  done  by  monitor¬ 
ing  the  error  rate  from  pilot  information  placed  in  message  headers. 

It  is  also  observed  from  Figure  3.2  that  the  improvement  for  DDE 
over  NDDE  increases  with  diversity  order.  This  is  because  as  the 
number  of  diversity  branches  increase,  so  does  the  excess  noise  due  to 
NT)DE. 

Figure  3.4  compares  the  use  of  the  sub-optimal  sample  mean  esti¬ 
mator  versus  the  Wiener-Hopf  estimator  of  Section  3.3.1.  The  com¬ 
parison  was  for  a  uniform  delay  power  profile,  DDE  and  =  10.  It 
is  seen  that  over  the  range  of  practical  error  rates  (P(e)<10“^)  the 
more  practical  sub-optimal  estimator  performs  as  well  as  the  Wiener- 
Hopf  estimator.  At  low  SNR  the  improvement  obtained  from  the 
Wiener-Hopf  estimator  is  only  observed  for  higher  orders  of  diversity. 
At  low  values  of  SNR,  to  achieve  the  same  performance  as  the  optimal 
Wiener-Hopf  estimator,  the  sample  mean  estimator  must  use  a  larger 
memory  length. 

Figure  3.5  illustrates  the  effect  of  increasing  the  memory  length  of 
the  estimator.  It  is  obsei-ved  that  for  a  channel  with  uniform  delay 
power  profile  as  M^,  — ►  oo  the  probability  of  error  approaches  the 
lower  bound,  due  to  perfectly  coherent  reception.  It  is  observed  that 
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Figure  3.2  Probability  of  error  for  DDE  and  NODE  -  Uniform  delav  power 
profile. 
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Figure  3.3  Probability  of  error  for  DDE  and  NODE  -  Exponential  delay  power 
profile. 
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Figure  3.5  Probability  of  error  for  increasing  lengths  Figure  3.6  Effect  of  diversity  order  on  memory  length 
of  estimator  memory  -  Uniform  delay  -  Gaussian  delay  power  profile 

power  profile. 


Figure  3.7  Ciomparison  of  approximation  and  exact 
expressions  for  probability  of  error  of  par¬ 
tially  coherent  receiver  -  Exponential  delay 
power  profile. 


Figure  3.8  Effect  of  decorrelated  estimates  on  proba¬ 
bility  of  error  -  Uniform  delay  power 
profile. 


for  Mjj  =  10  the  degradation  in  performance,  for  P(e)  <  10“^,  is  less 
than  0.5  dB.  Short  memory  lengths  are  preferred  to  present  run-away 
errors. 

Figure  3.5  and  Figure  3.7  illustrate  that  although  an  increase  in  ' 

order  of  diversity  improves  performance,  it  is  necessary  to  increase  the 
memory  length  simultaneously.  A  useful  interpretation  relating 
increasing  memory  length  and  diversity  order  is  obtained  as  follows. 

It  is  desirable  to  choose  the  memory  length  of  the  sample  mean  estima¬ 
tor  so  that  it  has  the  same  error  variance  as  the  Wiener-Hopf  estima-  j 

tor. 

Equating  (3.27)  and  (3.29)  gives: 

,  _ ^  (3,5S, 

Tp  4'.(kTp)  Mb  +  hN„  A^EpMb' 


where  is  the  memory  length  of  the  sample  mean  estimator  that  is 
required  to  obtain  the  same  error  variance  as  the  optimal  VViener-Hopf 
estimator,  with  memory  length  Now  (3.58)  can  be  rewritten  as: 


Mb  =  Mb  + 


hNo 

Tp  ^kT,) 


(3.59) 


Hence,  for  a  fixed  error  rate  (fixed  N^)  and  fixed  code  rate  (Tp),  Mj 
depends  inversely  on  the  average  path  strength,  ^c(kTp).  For  large 
orders  of  diversity,  the  inclusion  of  weak  paths  cause  Mb  to  get  very 
large.  This  increase  in  memory  length  should  not  exceed  the  channel 
coherence  time  (T^),  since  if  MbT  >  T^,  then  the  estimator  filter  will 
have  a  bandwidth  smaller  than  the  fade  rate  and  decorrelated  esti¬ 
mates  will  be  obtained  followed  by  performance  degradation.  This 
consideration  may  then  dictate  lower  orders  of  diveisity,  so  as  to  avoid 
including  weak  paths.  Figure  3.7  also  illustrates  the  tightness  of  the 
approximation  given  by  (3.46)  and  (3.48),  and  hence  its  utility  in 
evaluating  the  error  probability  of  a  partial  coherent  receiver  for  a 
non-uniform  profile. 


Figure  3.8  compares  the  upper  bound  using  the  high  SNR  approx¬ 
imation  and  decorrelated  estimates,  for  =  10“®  (~  0),  with  the 
non-coherent  receiver  performance.  The  agreement  between  the  two 
curves  indicate  that  when  the  channel  conditions  preclude  estimation, 
it  is  not  possible  to  do  better  than  non-coherent  detection. 


3.7  Receiver  Structures  for  Polyphase  Signaling 


The  received  signal,  in  the  symbol  interval  (L— 1)T  <  t  <  LT, 
under  the  two  hypotheses  (H,,  i  =  0,1)  is  given  by: 

H,  =  H(L) 

■■(t)  -  S  Ci(L)  A  p(t  -  kTp)  +  Z(t)  (3.C0) 

k-1 
Ott 

where  ^j(L)  =  i  =  0,  1.  p(t)  is  the  DSSS  code  and  Z(t)  is  the 

AWGN  process  with  psd  of  N^. 

The  receiver  structure  is  derived  as  in  Section  3.2.  The  receiver  struc¬ 
ture  that  decides  between  the  two  hypotheses  (Hj,  Hq)  in  the  signaling 
interval  (L— 1)T  <  t  <  LT  can  be  shown  to  be: 


Nc 

E 


NT  .  H(L)  -  H, 

--  A-r'o  '  ..v  R^[»-k(L)  C,*(L)]  >  0 

k_i  A-E  2ct£,,  +  N^,  H(L)  -  H. 


(3.61) 


where 


LT 


i-k(L)  =  /  rk(t)  — = 

(L-i)T  VEp 


p(ti  -  kTp)  dt  ;  tj  =  t  -  (L-l)T  (3.62) 


Ck(L)  and  are  defined  by  (3.11)  and  (3.12).  As  in  Section  3.2  the 
estimate  Ck(L)  is  based  on  the  memory  information 
^kL  =  [Vk(L-Mb)  •  •  •  Vk(L-l)]'^,  where  Vj^(j)  is  given  by: 


^kO)  =  dj  ^kO) 


(3.63) 


dj,  is  the  bit  decision  in  the  signaling  interval.  dj=l  when  Hj  is  true 
and  dj=0  when  is  true.  As  before,  the  total  memory  information  is 
given  by: 


II"  =  l^iL  m2L 


m 


n^l] 


(3.64) 


where  nriiL,  ni2L  '  ’  ’  are  independent  blocks  of  Gaussian  memory 
information,  each  block  being  used  to  estimate  one  path  coefficient. 

From  (3.72)  and  (3.73)  the  receiver  structure  can  be  interpreted  as 
an  adaptive  estimator-correlator  that  provides  diversity  through  Maxi¬ 
mal  Ratio  Combining,  (MRC).  The  structure  of  the  receiver  is  shown 
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in  Figure  3.9  and  Figure  3.10. 

Mavimal  Rntin  Combiner  (MRC) 

Assuming  that  the  sample  mean  estimator  is  used,  substituting  (3.29) 
and  (3.30)  in  (3.61)  yields: 


Nd 

S  Re 


k-i 


L-l 


j-L-Mb 


rk(L)  r,£*(j) 


H(L^Ho 

H(Lr“H, 


0 


(3.66) 


where  Nj  <  N^,  is  the  order  of  diversity  employed.  In  practice  Nj  is 
usually  chosen  to  be  less  than  the  total  number  of  resolved  paths,  for 
implementation  reasons. 

2.  Differential  Detection 

A  differential  detector  can  be  obtained  from  (3.66)  by  using  a  one 
bit  memory.  For  Mb=l  (3.66)  reduces  to: 

Nd  H(L)-Hi 

E  Re[rk(L)  rk*(L-l)]  >  0  (3.67) 

k-l  H(L)-Ho 


This  differential  detector  is  used  when  Differential  Phase  Shift  Keying 
(DPSK)  is  the  source  modulation. 

3.  M-ary  Detection 

For  binary  phase  coded  signaling  the  receiver  structure,  shown  in 
Figure  3.9,  can  be  considered  to  be  obtained  by  replacing  the  conven¬ 
tional  matched  filter/correlator  in  the  narrowband  system  with  struc¬ 
ture  of  (3.62).  Similarly,  for  M-ary  phase  coded  signaling,  the  receiver 
structure  can  be  obtained  from  its  narrowband  analog  by  replacing  the 
matched  filters  by  the  structure  given  in  Figure  3.10.  The  structure  of 
this  M-phase  detector  is  illustrated  in  Figure  3.11. 

3.8  Performance  of  Coherent  Receiver 

Proakis,  [19],  has  obtained  probability  of  error  expressions  for 
adaptive  reception  of  M-phase  signals  considering  diversity.  The 
derivation,  therein,  is  fairly  complex.  Further,  the  case  of  a  non- 
uniform  delay  power  profile  has  not  been  considered 

3.8.1  Uniform  Delay  Power  Profile 


Figure  3.9  Receiver  structure  for  BPSK  signaling. 


Figure  3.10  Estimator-correlator  structure  for  BPSK  si 


signaling. 


In  this  case  all  paths  are  similar  and  Tp  •  $(,(kTp;0)  =  for 
k  =  1,  2,  N(,.  Hence,  all  paths  have  the  same  average  SNR. 
Polyphase  signaling  is  used  and  the  receiver  structure  of  Figure  3.11  is 
used  to  demodulate  the  received  signal. 

The  basic  decision  variables  Ui(L),  and  Vj(L),  from  Figure  3.11 
and  (3.66),  can  be  written  as: 


Ui(L)  =  E  ReKAVE^ell'M  Cj(L)  +  Zt(L)}  C,*(L)1 

k-1 


(3.68a) 


Nd 


v,(L)  -  S  ReKAVE^e  ^  c„(L)  +  Zk(L)}  Ct‘(L)]  (3.68b) 

k-1 


where 


LT 


Zk(L)=  /  Z(t) p{t,  -  kTp)  dt  ;  t, -t-(L-l)T 

(L-l)T  Vhip 


{Zj5(L),  k  =  1,  2,  ...,  N(j}  are  a  set  of  i.i.d.  zero  mean  complex  Gaussian 
noise  random  variable  with  variance  and 

&,(L)  =  27r(i-l)/M,  i  »  1,  2,  ...,  M. 

The  probability  of  error  is  derived  by  first  determining  P(e  ll™)- 
If  the  memory  information,  I^,  is  given,  it  is  equivalent  to  knowing 
{Ci^(L),  k  =  1,  2,  •  •  •  Nd}.  In  order  to  determine  P(e  ll™),  the  statis¬ 
tics  of  the  decision  variables  Ui(L)  and  Vj(L)  conditioned  on  need  to 
be  known.  It  can  be  shown,  using  (3.68a)  and  (3.68b),  that  the  random 
variables  Ui(L)  and  Vj(L),  conditioned  on  the  signaling  phase  ^,(L)  and 
I^,  are  Gaussian  with  the  following  statistics: 


Nd 

E|u,(L)  I  Il".«,(L)1  -  aVe;cos«|(L)  E  |Ck(L)  P 

k-1 


Nd 

E|v,(L)  I  =  AVE;sine,{L)  E  |Ck(L)  P 

k-1 


Var|u,(L)  |  If',«,(L))  =  S  (A^c'S.  +  -^1  |Ck(L)  P 

k-1 
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Varlv,(L)  |  If  .«>,(L)|  =  £  1A%<Tl\  +  i^]  |Ct(L)  \- 

k-1  ^ 


The  problem  is  to  detect  one  of  M  equiprobable,  equal  energy 
polyphase  signals,  and  is  similar  to  that  described  in  Section  4-5  of 
[lO].  Hence,  the  conditional  probability  P(e  |l^)  can  be  obtained  from 
the  derivation  of  Section  4-5  replacing  the  corresponding  decision  vari¬ 
ables  by  {ui(L)  IIl',  ^i(L)}  and  {vj(L)  |Il",  ^i{L)}.  The  result  is: 


P(e|ir)  =  ^-|erfV^sin^ 


1 

-  /  exp(-y^)  erf  y  cot^  dy 


(3.69) 


where 


A=E„  S  |C,(L)P  2 


(3.70) 


S  lA«Ep  2<T£k  +  NJ  |C|,(L)  r- 

k-1 


Now  using  (3.29)  in  (3.70)  yields 


N„|l  + 


N.  . 

E  ICk(L)P 

k-1 


(3.71) 


where 


XL  =  E  |Ck(L)  1“  and  2<t2  =  [l  + 


Now,  from  (3.69)  and  (3.71)  the  average  probability  of  error  can  be 
written  as: 
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OO 


P(e)  -  /  P(e  1  !£>)  P(xi.)  dx;, 
0 


OO 


=  /P(e 


I  Xl)  P(xl)  dx'L 


(3.72) 


For  a  uniform  delay  power  profile  P{xl)  is  given  by  (3.61).  Substitut¬ 
ing  (3.61)  and  (3.69)  in  (3.72)  and  integrating  yields: 


Nd-i  2k 

< 

^  ~  Ml  Yj  ( }^  ) 

A 

k-o 

4 

V  / 

Vi 


k-0 


,  2 
1  - 


(3.73) 


where 


OO 


I(k,  M)  =  /  exp 
0 


y* 

Ml 


erf 


y  cot(^) 


dy 


and 


Ml  = 


sin^(^)  +  2cr“ 


,  2\l  =  2\l, 


Lk 


1  o  ^ 

^Lk 


From  (3.73)  it  is  observed  that  for  M=2  (BPSK)  a  closed  form  expres¬ 
sion  is  obtained  and  is  given  by: 

BPSK 


P(e) 


2 


■ 

✓  \ 

c 

Nd-i  2k 

1  -  Atf 

^  ~  Ml  Y  (  k  ) 

4 

k-O 

1  ^  J 

(3.74) 
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where 


= 


A=E„ 


A=Ep  +  Noll  + 


For  M  >  2  the  probability  of  error  evaluation  requires  one  numerical 
integration  due  to  I(k,  M). 

3.8.2  Non-Uniform  Delay  Power  Profile 


This  case  corresponds  to  dissimilar  paths  and  hence  the  average 
SNR  on  each  path  differs.  The  probability  of  error  is  expressed  as  in 

(3.72) ,  where  P(xl)  is  given  in  (3.-1C).  Hence,  using  (3.46)  and  (3.69).  in 

(3.72)  and  solving  the  integral  yields: 


P(e)  = 


M-1 

M 


1  — ^  tan  ^ 

71 


(3.75) 


where 


and 


f^k  = 


♦  o  TT 

Sin-—"  2Xi 
M 


o 

Lk 


o 

2^e-  = 


A-E 


P  L 


1  + 


Mr 


Now,  from  (3.75)  for  M=2  (BPSK)  the  probability  of  error  is  given 
by: 


P(e)  = 


2 


Nd 

E  ^Lk 

k-1 


A^Ep  2X^, 


Lk 


A^Ep  2Xl\  +  N,[l  -h  ^ 


Mr 


(3.76) 
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3.8.3  Approximations  and  Bounds 

For  M  >  2  it  is  necessary  to  perform  a  numerical  integration  in 
order  to  evaluate  the  error  probability  over  a  channel  with  a  uniform 
delay  power  profile.  It  was  determined  that  by  approximating 

I  !£")  ^  good  approximation  to  the  error  probability,  in 

closed  form,  -an  be  obtained. 

From  [10],  (3.69)  can  be  approximated  as: 


P(e  I  I™)  ~  erfc 


sin-^ 


M 


(3.77) 


Using  (3.77)  in  (3.72)  and  performing  the  appropriate  integrations,  for 
uniform  and  non-uniform  profiles,  the  following  expressions  for  proba¬ 
bility  of  bit  error  are  obtained: 

Uniform  Delay  Power  Profile 


P(e)  ~  1 


Nd-l  2k 

/^I  E  (k) 

k=l 


1  -  /'f 


M  >  2 


(3.78) 


M'here  Pj  is  defined  as  in  (3.73) 

Non-Uniform  Delay  Power  Profile 

Nd 

P(e)  -  E  '^Lk  A'k  >  2  (3.79) 

k=l 


where  is  defined  in  (3.75)  and  is  defined  in  (3.46). 

For  QPSK  (M=4)  and  Gray  encoding  of  source  bits,  it  is  known 
that  the  probability  of  error  is  given  in  [10]  as: 


Pqg(®) 


P(e) 

log2M 


(3.80) 


Substituting  (3.78)  for  P(e)  in  (3.80)  yields  (3.74),  indicating  that 
QPSK  with  Gray  encoding  of  source  bits  gives  the  same  error 
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probability  as  BPSK.  QPSK,  however,  has  twice  the  throughput  for 
the  same  error  rate. 

3.8.4  Differential  Detection 

The  probability  of  error  expressions  for  DPSK  and  differential 
detection  can  be  obtained  from  (3.75)  and  (3.77)  using  Mt>  ~  1- 

It  is  to  be  noted  that  all  the  probability  of  error  expressions  in 
this  section  hold  only  for  k  =  1,  2,  ...,  N^j.  This  is  rea¬ 

sonable  since  if  the  reverse  were  true  it  would  imply  that  the  estimated 
SNR  is  negative. 

3.9  Performance  Evaluation  for  Polyphase  Signaling 

The  performance  is  evaluated  with  regard  to  observing  the  effects 
of  delay  power  proffle,  diversity  order,  estimator  memory  size,  R.M.S 
delay  spread  and  signaling  alphabet  size  on  probability  of  error.  The 
performance  of  the  coherent  receiver  is  also  compared  with  that  of  a 
differential,  square-law  and  a  partially  coherent  receiver. 

The  performance  evaluations  are  conducted  for  a  quasi-static 
slowly  fading  Rayleigh  channel.  An  Uniform,  an  Exponential  and  an 
Irregular  delay  power  profile  were  used  in  evaluating  the  error  perfor¬ 
mance.  The  non-uniform  profiles  were  assumed  to  have  a  250  ns  delay 
spread  with  a  40  dB  dynamic  range  in  signal  levels.  The  Uniform 
profile  was  assumed  to  have  a  total  delay  spread  of  250  ns.  A  signal¬ 
ing  rate  of  32  kbps  and  a  code  length  of  255  chips  was  chosen.  It  is 
assumed,  as  before,  that  perfect  bit  and  chip  synchronization  exists. 
The  results  of  the  performance  are  illustrated  in  Figures  3.12  -  3.19 

The  following  conclusions  were  drawn  from  the  performance  evalua¬ 
tions: 

1.  The  partially  coherent  receiver  has  a  BER  that  is  bounded  on 
the  lower  end  by  the  BER  of  the  coherent  Rake  receiver  and  on 
the  upper  end  by  the  BER  of  a  Square-Law  Rake  receiver. 

2.  When  the  fade  rate  is  larger  than  the  bandwidth  of  the  estimator 
filter,  decorrelated  estimates  occur  and  the  BER  performance  of 
the  partially  coherent  receiver  approaches  the  upper  bound. 

3.  DDE  produces  a  lower  BER  than  NDDE  due  to  a  3  dB  gain  in 
phase  estimation  SNR,  under  the  ideal  assumption  that  no  bit 
decision  errors  are  made.  To  prevent  propagation  of  bit  errors  in 
memory,  it  is  suggested  that  short  memory  lengths  (M^,  <  10)  or 
exponential  weights  be  employed. 

4.  The  sample  mean  estimator  performs  as  well  as  the  Wiener-Hopf 
estimator  when  its  memory  length  is  chosen  by  (3.59). 
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Figure  3.12  Effect  of  delay  power  profile  on  probability 
of  error  for  coherent  receiver  (N4=3). 


Figure  3.13  Effect  of  delay  power  profile  on  probability 
of  error  for  coherent  receiver  (Nd=6). 


Figure  3.14  Probability  of  error  for  polyphase  signal- 

(Nj=2)  -  Exponential  delay  power 
profile. 


Figure  3.1.5  Probability  of  error  for  polyphase  signal¬ 
ing-  (Nj=4)  -  Exponential  delay  power 
profile. 


Pr-At>AitUtf  ot  t/rmr 


10* 


Logeno 


e  " 

-  :k 

•  ra  a 

J-w«9a  o<c»>o  <>M  Mr  M  (oei 


Figure  3.16  Probability  of  error  for  polyphase  signal¬ 
ing,  (Nd=8)  -  Exponential  delay  power 
profile. 


yf. 


Figure  3.18  Comparison  of  probability  of  error  for 
different  receivers  (Na=2)  -  Exponential 
delay  power  profile. 


Figure  3.17  Probability  of  error  v-ersus  normalized 
delay  spread  -  Gaussian  delay  power 
profile. 


Figure  3.19  Comparison  of  probability  of  error  for 
different  receivers  (Nj=4)  -  Exponential 
delay  power  profile. 
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5.  Large  orders  of  diversity  result  in  inclusion  of  weak  paths  and 
hence  large  memory  lengths  (3.59).  This  results  in  decorrelated 
estimates  and  an  increase  in  BER.  Hence,  lower  orders  of  diver¬ 
sity  are  recommended. 

6.  For  a  quasi-static  channel  coherent  reception  using  MRC  and 
DDE  is  feasible. 

7.  The  coherent  receiver  like  the  ONMD  receiver  is  sensitive  to  the 
shape  of  the  delay  power  profile.  This  emphasizes  the  need  for 
measuring  indoor  channel  profiles  to  determine  the  appropriate 
transmission  power. 

8.  The  coherent  receiver  is  sensitive  to  R.M.S  delay  spread  variation 
and  exhibits  a  greater  degradation  with  increasing  orders  of 
diversity.  This  also  suggests  lower  order  of  diversity.  The 
receiver  is  more  sensitive  to  delay  spread  than  the  shape  of  the 
delay  power  profile. 

9.  Considerable  diversity  gain  is  achieved  in  using  higher  orders'  of 
diversity.  This  is  offset  by  the  sensitivity  to  R.M.S  delay  spread 
and  increase  in  memory  length.  An  encouraging  alternative  to 
achieve  higher  ordere  of  diversity  is  to  use  a  low  order  of  mul¬ 
tipath  diversity  and  employ  multiple  antennas. 

10.  The  use  of  a  higher  signaling  alphabet  results  in  a  severe  perfor¬ 
mance  degradation  for  M  >  4.  BPSK  and  QPSK  with  Gray 
encoding  or  source  bit  perform  equally  well.  QPSK,  however, 
offers  twice  the  data  rate  within  the  same  bandwidth. 

11.  For  a  quasi-static  channel  the  use  of  phase  coded  signaling 
(BPSK)  and  MRC  results  in  3  dB  improved  performance  over 
differential  detection  using  DPSK  signaling. 

12.  Orthogonal  Signaling  using  partially  coherent  receptic>n  performs 
as  well  as  DPSK  using  differential  detection.  However,  the  sim¬ 
plicity  of  the  differential  detector  favors  its  use  in  practice. 
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Probability  of  error  for  QPSK  and  DPSK  in  a  CDMA  system. 


Effect  of  increase  in  code  length  on  probability  of  error  for  Orthog 
onal  Signaling  (O.S.)  relative  to  BPSK. 
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Figure  4.7  Simulation  results  for  probability  of  error 
for  BPSK  considering  perfect  estimation 
(Na=2) 
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Figure  4.8  Simulation  results  for  probability  of  error 
for  BPSK  considering  estimation  errors 
(Nd=2). 
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Figure  4.9 


Simulation  results  for  probability  of  error 
for  BPSK  considering  estimation  errors 
(Nd=4). 
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Figure  4.10  Simulation  results  indicating  effect  of  imperfect  estimates. 
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Figure  4.11  Simulation  results  indicating  effect  of  interpath  interference. 
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CHAPTER  4 

COMPARISON  OF  MULTIPATH  DIVERSITY 
RECEIVERS  IN  A  CDMA  SYSTEM 


4.1  Introduction 

It  is  the  purpose  of  this  chapter  to  study  the  performance  of  three 
multipath  diversity  receivers  in  a  multiple-user  environment  (CDMA) 
taking  into  account  ISI,  interuser  interference  and  non-ideal  code  corre¬ 
lation  effects.  In  evaluating  the  probability  of  error  for  the  different 
receiver  structures,  the  following  CDMA  system  is  considered.  The 
channel  is  an  indoor  wireless  communication  channel.  Each  user  can 
transmit  speech  at  a  rate  of  32  kbps.  Data  will  be  ’  l  ansmitted  in 
packets.  It  is  assumed  that  the  user  terminals  communicate  with  the 
base  station  in  a  star  network  architecture,  [7],  [8].  The  users  com¬ 
municate  to  the  base  station  using  asynchronous  CDMA. 

The  central  station  is  assumed  to  have  a  bank  of  receivers,  one 
for  each  unique  user  code.  These  receivers  are  multipath  diversity 
receivers.  Each  user  terminal  has  a  similar  receiver  structure.  To 
avoid  the  near-far  problem  average  power  control  is  assumed. 

4.2  Performance  Analysis 

4.2.1  Coherent  Multipath  Diversity  Receiver 

Considering  BPSK  the  band-pass  transmitted  signal  for  the  k^^ 
user  in  the  CDMA  system  is  represented  as: 

Sk(t)  =  Re[A  Uk(t)  (4.1) 


where  f^,  is  the  carrier  frequency,  A  is  the  signal  amplitude  and 

Uk(t)  =  bk(t)  Pk(t)  (4.2) 

where  bk(t)  is  the  data  bit  stream  of  the  k*'^  user  and  Pk(t)  is  the 
corresponding  DSSS  code. 

bkW-  S  b.k  ET(t  -  iT)  ; 


(4.3) 
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OO 

Pk(0=  S  Pk[j]  gTp(t  -  >Tp)  ;  {p,,[j]}  =  tl  (4.4) 

J  — OO 


where  g^{t)  =  1  0  <  t  <  r  and  zero  otherwise.  It  is  assumed  that  the 

DSSS  codes  are  periodic  with  period  Np  and  that  T  =  Tp  Np.  The 
energy  in  each  period  of  the  code  sequence  is  denoted  by  Ep. 

Consider  K  asynchronous  users  of  the  communication  channel. 
Assuming  the  channel  model  of  [12]  the  received  signal  in  the  interval, 
(L— 1)T  <  t  <  LT,  can  be  written  as: 

Kt)  =  E  E  A  C^{L)  Uk(t  -  tk  -  mTp)  +  T]{i)  (4.5) 

k— 1  m— 1 

where  ??(t)  is  an  additive  white  zero  mean  complex  Gaussian  noise  pro¬ 
cess,  at  the  receiver  input,  with  autocorrelation  function 
<^„(t,  r)  =  Nq  6(t  —  r).  tjj  is  a  uniform  random  variable,  in 
[0,  T  —  T^,],  representing  the  time  of  transmission  for  each  user.  It  is 
assumed  that  k=l  corresponds  to  the  reference  user  and  that  tj  =  0. 
C^(L)  is  the  m*'^  channel  coefficient  for  the  user  in  the  interval 
(L— 1)T  <  t  <  LT.  The  C,|^(L)  are  uncorrelated,  zero  mean  complex 
Gaussian  random  variables  with  variance 

E|  |C^(L)  PI  -  2Qk=„(L)  =  Tp  •J>‘(mTp;0)  (4.6) 


where,  <l>p’'’(T;t)  is  the  channel  delay  power  profile  for  the  k*'^  user. 

The  receiver  structure  is  given  by  (3.66)  and  can  be  rewritten  for 
the  reference  user  (k=l)  as: 


S  Re|r,>(L)  C,'  (L)| 

J-1 


H(L)  -  H, 
H(L)  -  Ho 


0 


(4.7) 


where 


LT 


rm(L)  =  /  r(t) 

(L-l)T  VEp 


Pi(ti-mTp)dt  ;  t,  =  t  -  (L-l)T  (4.8) 


and  C^(L)  is  the  MMSE  estimate  obtained  as  follows: 
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1 

aVE^  Mb 


E  V'-mO) 

j-L-Mb 


(4.9) 


where  bj  is  the  bit  decision  in  the  j’'*'  signaling  interval  for  the  refei- 
ence  user. 

Substituting  (4.5)  in  (4.8)  and  recalling  that  the  codes  are 
periodic,  the  output  of  the  correlator  on  branch  of  the  receiver  can 
be  written  as: 


A  T 

V^p  m=l  0 

A  K  Nd  , 

+  T-T^  EE/  Cm(L)  Pk(t  -  ^mk)  bk(t  -  r^k)  Pl(t  “  jTp)  dt 

v'  ^p  k=2  m=l  0 
LT 

VEp  (L-l)T 


where  =  t^.  -1-  mTp  which  is  a  uniform  random  variable  in  the 
interval  [0,T].  Recall  rjj,j  =  mTp.  Substituting  (4.3)  in  (4.10)  and 
using  the  method  outlined  in  [20],  [6]  (4.10)  can  be  expressed  as: 

rj'{L)  =  A  Ve^  bi!  C,>(L)  +  X,(L)  +  Y,(L)  +  J)/(L)  (4.11) 

where 


X,(L)  = 


N, 


VeT 


s  Cj.(L)[b;  R„(4,)  +  b^.,  R„(r>  )] 


p  m— 1 

nifij 


Y,(L)  - 


N, 


VeT 


E  Cji(L)(b,;  RK,(rk,)  +  b,;.,  Rk,(rk,)l 


p  m«l 
m^l 


Rki(^)  =  /  Pi(t)  Pk(t  -  r)  dt  k  =  1,  2,  •  •  •  K 


r 
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Rki(^)  =  I  Pi(t)  Pk(t-^)  dt  k  =  1,  2, 


=  ^-k  +  mTp  -  j  Tp  ;  k  =  1,  2,  •  •  •  K  ;  =  0 


where  Rki(^)  and  RkiCT")  are  continuous  partial  correlation  functions. 


Let  Zj(L)  =  Xj(L)  +  Yj(L)  +  r}^\^  represent  the  total  interference  on  the 
j*'*’  branch  of  the  reference  user.  Xj(L)  represents  the  interpath 
interference  arising  due  to  the  non-ideal  autocorrelation  of  the  DSSS 
codes.  Yj(L)  represents  the  interuser  interference  arising  due  to  the 
non-ideal  cross-correlation  of  the  DSSS  codes. 

Substituting  (4.11)  in  (4.7)  the  basic  decision  variable  becomes: 


Ui(L)  =  E  Re'  \a  bL^  C/(L)  +  Zj 


(L)]c,'’(L) 


(4,12) 


The  following  assumptions  are  now  made 

(i)  The  delay  power  profile  is  uniform. 

(ii)  Xj(L)  is  small  in  comparison  to  Yj(L).  This  is  particularly  true 
for  low  orders  of  diversity  and  a  large  number  of  active  users. 

(iii)  Zj(L)  =  Vj(L)  +  is  a  Gaussian  random  variable.  This  is  a 
logical  assumption  based  on  the  central  limit  theorem.  It  can 
also  be  shown  that  Zj(L)  is  zero  mean  and  has  a  variance 
which  is  computed  below. 

(iv)  E[Zj(i)  Zj*(k)]  =  0  when  i?^k. 

(v)  E[Zj(i)  Zni*(i)]  =  0  when  j  ^  m. 

Under  these  assumptions  the  decision  variable  given  by  (4.12)  is 
similar  to  that  given  by  (3.68a).  Hence,  the  corresponding  probability 
of  error  will  be  the  same,  except  that  Nq  is  replaced  by  which  is  the 
variance  of  the  new  noise  variable  Zj(L). 

Assuming  that  the  interpath  interference  is  much  smaller  than  the 
interuser  interferene,  it  can  be  shown  that: 


=  Et 


,  K  Nd 

r  S  E  2a|?„(L)  +  N„ 

T  k-2  m-1 


(4.13) 
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where  E^,  =  A'T  is  the  transmitted  energy  per  symbol  period.  For  a 
uniform  profile,  2Qirjjj(L)  =  .  m  =  1,  2,  and 


= 


(E14) 


4 


Hence,  for  a  non-uniform  delay  power  profile,  the  probability  of  error 
is  given  by  (3.77)  with  replaced  by  which  is  defined  in  (4.13). 
For  a  uniform  delay  power  profile  the  probability  of  error  is  given  by 
(3.75)  with  N,.  replaced  by  a^,  which  is  defined  in  (4.14).  For  M  >  2, 
tne  probability  of  error  expressions  given  by  (3.73)  and  (3.76)  can  be 
used,  appropriately,  with  replaced  by 

4.2.2  Differential  Multipath  Diversity  Receiver 

The  receiver  structure  for  the  user  may  be  expressed  as; 

.Nd  ,  ,  H(L)  -  H, 

s  R4'k(L)  „  0  (4.15) 

m=l  "I*-)  “  Ho 


i 


n 


i 


i 


This  receiver  structure  can  be  obtained  from  (4.7)  and  (4.9)  using 
Mt,  =  1.  Hence,  the  probability  of  error  for  the  differential  receiver  can 
be  obtained  from  the  probability  of  error  expressions  of  the  coherent 
receiver  by  using  M  =  2  and  Mb  =  1. 

4.2.3  Non-Coherent  Multipath  Diversity  Receiver 

Each  user  employs  a  pair  of  orthogonal  codes  to  transmit  binar>' 
information.  Let  Pt^i(t)  and  Pko(t)  denote  the  pair  of  codes  used  by  the 
user.  The  received  signal  can  be  represented  as: 

r(t)  =  E  E  A  C^(L)  §  a,*"  gT(t  -  “  >T)  Pki(t  -  T-mk) 

k— 1  m—1  I— — oo 


I 


I 


I 


+  (1  -  a,'')  gT(t  -  T-mk  -  >T)  Pko(t  -  7'ink)  +  v{t) 


(4.16) 


I 


Pk,  =  E  Pkilj]  Sxjt  -  j  T  )  ;  {Pk,!j]}  =  +  1  :  i  =0.  1 


(-1.17) 


gr(t)  =  1  0  <  t  <  T  and  zero  otherwise.  A  is  the  signal  amplitude, 

a,^  indicates  which  hypothesis  is  true  in  the  i*^^  signaling  interval. 
aj’^MD  when  H(i)=Ho  and  a,*'=l  when  H(i)=Hi.  tj.  is  a  uniform  random 
variable,  in  [0,T— T^j,  representing  the  random  time  of  transmission  of 
each  of  the  K  asynchronous  useis.  r?(t)  is  a  zero  mean  complex  Gaus¬ 
sian  random  process  with  autocorrelation  function  S(t~r). 

Assuming  a  uniform  delay  po%ver  profile,  the  receiver  structure  for 
the  user  is  expressed  as: 


X.t  HsL)  =  H,  Nd 

E  „  E  hfn.(L)|- 


(4.18) 


where 


fkm(L)  =  /  r(t)  pjti  -  mTp)  dt 

(L-l)T 


i  =  0.  1 


(4.19) 


and  tj  =  t  —  (L— l)T 

Assume  that  in  the  L'*'  signaling  interval,  the  reference  user 
transmits  Pn(t).  Now,  substituting  (4.16)  in  (4.19)  it  can  be  shown 
that  the  output  of  the  branch  correlator  of  the  reference  user's 
receiver  is: 


(L)  =  aV^  C,'(L)  +  X„(L)  +  Y,,(L)  +  »,>,  (L) 


(4.20) 


^  ^io(O) 

(L)  =  aVe;  c/(L)d-X,(L)-bY,/L)-f  7/,'/L) 


(4.21) 


where,  for  i  =  0.  1 


X,(L)  =  E  C/(L)  ail  R,V(r^)  +  (1  -  a^)  R,)/(r^., 
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+  “L-l  +  (1  -  ai!.,)  R„,(ri) 

~  ^  ^  [^L  (^mj)  +  (1  ~  ^l)  ^oV(’’mj) 

V^p  k-2  m-1  '• 

+  RiV«)  +  (1  -  ^Li)  RoV  (^i) 

T 

Rh  ^(7-)  =  /  Pkl(t-7-)  Pi:(tj  dt  ;  1  =  0,  1 

T 

T 

Rh  =  /  Pki(t-')  Pi,(t)  fit  ;  I  =  0,  1 
0 

=  tk  +  mTp  -  jTp  ;  k  =  1,  2,  •  •  •  K  ;  tj  =  0 


LT 

^ij(L)  =  -  /  ty(t)  Pi,(ti  —  jTp)  dt 

VEp  (L-i)T  I 

^io(O)  =  E  PiiPl  PioDl 

j-o  ! 


^io(O)  is  the  periodic  cross-correlation  of  the  DSSS  codes.  Also,  a,''  =  1 
ii*  Pki(t)  is  sent  in  the  i*'^  signaling  interval  and  a,*'  =  0  if  Pko(t)  is  sent 
in  the  i*^^  signaling  interval;  the  choice  determined  as  a  Bernoulli  event. 
R,,  (r)  and  Ril'n’')  are  continuous-time  partial  correlation  functions 

^ij(E)?  i  =  0,  1  represent  two  zero  mean  complex  Gaussian  random 
variables  with  variance  N^.  The  correlation  between  these  noise  vari¬ 
ables  can  be  shown  to  be: 


I 
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E[j//i(L)  =  X. 


(4.22) 


For  each  receiver  the  total  interference  in  the  t>  i.  arms  (i  =  0.  1)  are 
defined  as: 


Z,(L)  =  X,j(L)  +  Y„(L)  +  >J,|(L)  .  i  =  0.  1  (4.23) 


where  X,j(L)  is  the  interpath  interference,  V,j(L)  is  lue  interuser 
interference  and  t?,J(L)  is  the  AWGX  interference. 

Now.  as  in  Section  4.2.1  the  following  assumptions  are  made: 

(i)  The  delay  power  profile  is  uniform. 

(ii)  The  interpath  interference  is  negligible  in  compaidson  to  the 
interuser  interference. 

(iii)  Z|j(L)  =  X,j(L)  +  t/|j(L)  is  modeled  as  zero  mean  complex  Gaus¬ 
sian  noise,  with  variance  (T*,  from  the  central  limit  theorem. 

(iv)  E[Z,j(L)  Z,j*(k)]  =  0  when  L5>^k 

(v)  E[Z,n,(k)  Z,n*(k)]  *  0  when  m^n 


(vi) 


^/o(0) 


.X 


1 


Using  the  above  assumptions  it  can  be  shown  that  r/j(L)  and 
rj°(L)  are  zero  mean  correlated  complex  Gaussian  random  variables 
with  variances  <7f  and  cr^.  The  variance  of  each  of  these  random  vari¬ 
ables  can  be  shown  to  be: 


O 


1  -b  Nd(K-l) 


+  No 


(4.24) 


O 


Wo(Q)F 


+  Nd(K  -  1)  : 


3X% 


+  X'o 


(4.25) 


where  =  A^T  is  the  transmitted  energy  per  symbol  period. 
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Further,  it  can  be  shown  that  the  correlation  between  the  vari¬ 
ables  ri?jn(L)  and  ri?jj,(L)  is  given  by: 


>‘lo  =  E|,-^„{L)  rg„*(L)) 


^lo(O) 


N. 


(4.26) 


To  compute  the  probability  of  error,  the  decision  variable  from  (4.18) 
is  expressed  as: 


r  = 


E 

1=1 


',(L)|--  |r?,(L)P] 


(4.27) 


and  the  probability  of  error  is  given  bj’.* 


P(e)  =  P(r  <  0) 


(4.28) 


Using  equation  (4.21)  of  [18]  and  making  appropriate  changes  based  on 
the  statistics  of  i‘mk(L)  and  r^i((L)  it  can  be  shown  that: 


P(e)  = 


1  + 


V2 

V, 


N,j-i  2  Nrf— 1 

PXFT  S  ( 

m=0 


m 


V, 


(4.29) 


where, 

u  =  af  o-Q-  -  l/i/o  P 


\V  = 


o 


2u 


V,  =  (W2  -f  l-)''2  _  W 
u 

V2  =  (W2  -t-  -L)1/2  +  w 


4.3  Receiver  Performance  Comparison 


In  this  section  a  comparison  is  made  of  the  BER  performance  of 
three  receivers,  under  the  following  conditions:  uniform  delay  power 
profile  with  Tjjj  =  400  ns,  Rj  =  32  kbps,  Np  =  255,  Efj/Nj,  =  — 10  dB, 
K  =  15,  N(,  =  4  and  Mj,  =  10. 

Figure  4.1  illustrates  the  effect  of  a  non-ideal  cross-correlation  for 
Gold  codes  used  in  binary  orthogonal  signaling.  It  is  seen  that  the 
code  cross-correlation  raises  the  irreducible  BER  due  to  a  correlation 
between  the  diversity  branches.  Figure  4.2  illustrates  that  for  Kasami 
codes  the  degradation  from  ideal  performance  is  small.  This  is  because 
these  Kasami  codes  are  optimal  [18],  as  they  satisfy  Welch’s  lower 
bound  for  the  peak  cross-correlation  value. 

Figure  4. 3-4. 5  illustrate  a  comparison  of  three  modulation  and 
detection  schemes.  It  is  observed  that  interuser  interference  produces 
an  irreducible  BER  or  error  floor.  From  the  error  floors  it  can  be  seen 
that  BPSK-coherent  detection  offers  almost  a  3dB  advantage  over 
DPSK-differential  detection.  DPSK-differential  detection  offers  almost 
a  4  dB  advantage  over  orthogonal  signaling  (O.S.)-  non-coherent  detec¬ 
tion.  It  is  also  seen  that  for  higher  orders  of  diversity,  =  8,  there  is 
a  significant  performance  improvement.  Since  =  4,  this  order  of 
diversity  can  be  achieved  through  antenna  diversity  by  employing  two 
antennas.  Figure  4.5  illustrates  an  interesting  result,  where  it  is  seen 
that  DPSK-differential  detection  produces  a  lower  irreducible  BER 
than  QPSK-coherent  detection,  indicating  that  DPSK-differential 
detection  performs  better  when  the  interuser  interference  becomes  the 
limiting  noise  source. 

Non-coherent  detection,  although  simple  to  implement,  offers  very 
poor  performance.  To  obtain  better  performance,  for  a  given  order  of 
diversity,  there  ai-e  two  alternatives.  The  first  is  to  use  a  longer  DSSS 
code,  results  for  which  are  shown  in  Figure  4.6.  It  is  seen  that  in  order 
to  obtain  a  performance  on  the  order  of  BPSK-coherent  detection,  it  is 
necessary  to  use  a  code  that  is  four  times  as  long.  However,  if  this 
performance  improvement  is  to  be  achieved  for  the  same  order  of 
diversity  the  data  rate  must  be  divided  by  4.  In  other  words,  for  a 
given  bandwidth  and  delay  spread,  a  performance  gain  can  be  achieved 
using  a  longer  code  but  the  data  rate  must  be  appropriately  reduced. 

4.4  Simulation  Study 

In  evaluating  the  performance  of  the  receivers  for  a  multi-user 
system  a  number  of  assumptions  were  made  for  tractable  analysis. 
These  assumptions  could  result  in  highly  biased  results. 

In  practice,  the  interpath  interference  and  the  interuser  interfer¬ 
ence  produce  some  correlation  between  the  diversity  branches  which 


-  64  - 


negates  the  effect  of  diversity  combining,  [2],  [9].  Further,  the  effect  of 
the  total  interference  may  be  such  that  the  simple  sample  mean  esti¬ 
mator  may  not  be  adequate  for  coherent  detection. 

In  order  to  observe  the  severity  of  these  effects  and  verify  the  ear¬ 
lier  assumptions  a  small  CDMA  system  was  simulated.  The  simula¬ 
tion  study  was  conducted  for  a  15  user  system  using  255  chip  long 
Kasami  sequences.  For  this  system  the  following  parameters  were 
assumed: 

(i)  A  uniform  delay  power  profile  with  a  multipath  spread  of  400 
ns.  This  allows  a  maximum  diversity  order  of  Nj  =  4  fur 
Kasami  codes  of  length  255  chips  and  a  data  rate  of  32  kbps. 

(ii)  Average  power  control  is  assumed. 

(iii)  Each  user  employs  BPSK  signaling  at  32  kbps  and  coherent 
detection  using  MRC. 

(iv)  The  sample  mean  estimator  of  (4.40)  with  M^  =  10  was  used  in 
all  receivers.  Since  the  probability  of  error  of  interest  is  around 
10“'*  the  effect  of  bit  decision  errors  on  DDE  were  neglected. 

(v)  The  Kasami  sequences  used  are  Auto-Optimal  Least  Side  Lobe 
Energy  (AO/LSE)  sequences,  [2l].  AO/LSE  codes  produce  very 
low'-interference  when  used  in  CDMA  applications.  These  codes 
were  generated  using  the  initial  loadings  given  in  [21]. 

The  Monte  Carlo  technique  was  used  to  simulate  the  system.  A 
sample  size  of  was  used  to  estimate  error  probabilities  of  the 

order  10“^.  From  [3],  it  can  be  shown  that  this  sample  size  produces  a 
99%  confidence  interval  of  [1.29  P(e),  0.77  P(e)];  where  P(e)  denotes 
the  estimated  error  probability. 

The  results  of  the  simulation  are  presented  in  Figures  4.7-4.11. 
The  simulation  runs  resulted  in  a  normalized  error  of  less  than  11%. 
The  normalized  error  variance  was  computed  to  be  less  than  0.03%. 
Based  on  the  limited  simulation  results  the  follow'ing  conclusions  are 
made: 

(i)  The  interuser  interference  can  be  modeled  as  additive  Gaussian 
noise  when  the  product  of  the  number  of  active  users  and 
resolved  paths  is  large. 

(ii)  The  interpath  interference  is  small  in  comparison  to  the  interuser 
interference  for  small  delay  spreads.  Hence,  it  can  be  neglected 
in  order  to  simplify  the  analysis. 

(iii)  The  sample  mean  estimator  performs  adequately  in  the  presence 
of  interuser  interference. 
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Figure  4.1 


Bounds  on  irreducible  error  for  Orthogonal 
Signaling  (O.S.)  using  Gold  codes. 


Figure  4.2 


Bounds  on  irreducible  error  for  Orthogonal 
Signaling  (O.S.)  using  Kasami  codes. 
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Figure  4.3  Probability  of  error  for  BPSK  and  DPSK 
in  a  CDMA  system. 
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Figure  4.4  Probability  of  error  for  DPSK  and 
Orthogonal  signaling  (O.S.)  in  a  CD\L4 
system. 
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4.5  Conclusions 

Three  multipath  diversity  receivers  using  coherent,  differential 
and  non-coherent  detection,  respectively,  were  analyzed  based  on  the 
assumptions  that  there  is  average  power  control  and  interuser  interfer¬ 
ence  is  Gaussian.  The  performance  results  indicated  that  BPSK  using 
coherent  detection  performs  the  best,  followed  by  DPSK  using 
differential  detection.  The  use  of  orthogonal  signaling  and  non-coherent 
detection  requires  the  use  of  long  codes,  which  in  turn  stipulates  lower 
data  rates  for  a  given  bandwidth  and  delay  spread.  It  is  seen  that 
antenna  diversity  in  conjunction  with  low  orders  of  multipath  diversity 
produce  significant  BER  performance  gain.  The  simplifying  assumptions 
involved  in  the  performance  analysis  were  verified  through  simulation. 
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